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In  this  book  we  present  features  of  calcu¬ 
lating  the  thermodynamic  functions  of  chemically 
reacting  systems  used  in  cycle  calculation.  We 
calculate  the  thermodynamic  diagrams  of  chemically 
reacting  systems.  We  examine  possible  schemes 
for  thermodynamic  cycles  in  chemically  reacting 
systems.  We  give  a  method  of  calculating  the 
regenerator  effectiveness.  We  examine  the 
influence  of  the  basic  parameters  on  the 
efficiency  of  power-plant  cycles.  We  present 
the  thermodynamic  optimization  of  gas  and  gas- 
liquid  cycles. 

The  book  contains  6  tables,  35  figures,  and 
38  references. 

The  book  is  designed  for  scientific  workers 
and  engineers  studying  questions  of  the  thermo¬ 
dynamic  analysis  and  selection  of  working  media 
for  power  plants . 
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INTRODUCTION 


The  problem  of  a  working  medium/coolant  at  the  present  time  i3 
one  of  the  pressing  problems  of  both  thermal  and  atomic  power  engi¬ 
neering.  New  working  media  and  coolants  are  being  sought  both  for 
temperatures  and  pressures  which  have  already  been  realized  in  power 
engineering  as  well  as  for  future  plants.  This  work  is  devoted  to 
studies  of  the  thermodynamic  cycles  using  chemically  reacting  working 
media. 

The  development  of  fixed  and  transportable  power  engineering, 
including  nuclear,  is  based  mainly  on  the  use  of  steam  as  the  working 
medium.  'Recent  technical  and  economic  studies  have  shown  that  in 
thermodynamic  cycles  using  steam,  the  economically  optimum  parameters 
and  limiting  efficiencies  have  been  achieved,  [1-4].  In  addition,  the 
physicochemical  properties  of  steam  ere  such  that  there  are  consid¬ 
erable  limitations  on  the  unit  power  of  single-shaft  turbine  units 

[4]. 

In  this  regard,  at  the  present  time  there  is  an  urgent  need  for 
finding  new  coolants  and  working  media.'  This  problem  is  particularly 
pressing  for  atomic  power  engineering  with  nuclear  reactors  using 
fast  neutrons  [5-7]. 

One  of  the  ways  for  solving  the  given  problem  might  be  the  use, 
as  working  media  and  coolants,  of  chemically  reacting  mixtures  in 
which  there  are  reversible  chemical  reactions  accompanied  by  the 
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heat  effects  of  chemical  reactions  and  a  change  in  the  number  of 
moles  [8-12].  An  increase  in  the  number  of  moles  during  heating  and 
a  decrease  during  cooling  result  in  increased  efficiency  of  the 
thermodynamic  cycle,  while  the  presence  of  the  heat  of  chemical 
reactions  assures  high  thermophysical  properties  of  the  mixture  due 
to  transfer  of  a  significant  quantity  of  heat  by  concentration 
diffusion  [8]. 

The  class  of  chemically  reacting  mixtures  is  large;  certain  of 
the  mixtures  examined  by  us  are  given  in  Table  1. 

Of  the  large  class  of  chemically  reacting  mixtures  in  which 
reactions  occur  with  an  increase  in  the  number  of  moles  and  the 
absorption  of  heat  during  dissociation,  the  chemically  reacting 
mixture  Nj04*i2\Oa**2NO+Ot  has  been  studied  most  completely. 

To  study  the  efficiency  of  thermodynamic  cycles  we  must  know  the 
change  in  thermodynamic  properties  over  a  broad  range  of  temperatures 
and  pressures.  There  have  been  practically  no  methods  developed  for 
calculating  the  thermodynamic  properties  of  chemically  reacting 
mixtures  [8].  The  occurrence  of  chemical  reactions,  however,  imposes 
its  own  specific  nature  on  the  calculation  of  thermodynamic  properties; 
therefore,  the  obtained  thermodynamic  dependences,  from  the  assumption 
of  constant  composition,  cannot  always  be  used  when  calculating  the 
properties  of  chemically  reacting  mixtures.  The  specific  nature  of 
the  behavior  of  the  thermodynamic  properties  imposes  its  features  in 
the  calculation  of  the  parameters  of  thermodynamic  cycles  using 
chemically  reacting  mixtures  as  the  working  medium. 

!  •  -  ■ 

In  this  book  we- present  the  methodology  for  calculating  the 
thermodynamic  properties  of  chemically  reacting  mixtures  and  thermo¬ 
dynamic  cycles  using  chemically  reacting  mixtures  as  the  working 
media.” 

The  book  was  written  using  data  from  research  at  the  Institute 
of  Nuclear  Power  Engineering  of  the  Academy  of  Sciences  of  the 
Belorussian  SSR  during  the  period  from  1964  thru  1969. 

1  I 
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Table  1.  Chemically  reacting  mixtures. 


(1) 

li 

(3) 

TcMaqvTyp* 

til1 

li 

(w)  ■  liquid 

J.o 

300-450 

13.7 

2NO,^2NO-fO, 

1.8 

400-1300 

27.0 

3.0 

300-1300 

13.7/27 

AIA^AJCJ, 

2.0 

475-1300 

29.8 

AljBr.-JAlBr, 

2.0 

600-1600 

30,0 

AM.-2AU, 

2.0 

500-1600 

26.4 

2NOBr^INO+Br, 

I.S 

300-600 

- 

2N0CUTJN0+O, 

l.s 

300-1200 

- 

Hed,rHC+a,5H|+ja 

s.o 

1300-3000 

- 

3.0 

900-1500 

- 

urfl 

2.0 

790-1700 

- 

CH.+2HiO-OOb-MH, 

2,0 

700-1500 

- 

.  CH,fO),-X»+2H, 

2.0 

600-1600 

— 

c.h.-ca+h, 

2.0 

700-1500 

- 

A1,Q,  f-  4AK*)iviAICl 

S.o 

930-1300 

263,8 

Al,Br,+  4Al<*)28AIBr 

6,0 

930-1600 

282,4 

AI,I,+4AK*)r#AII 

6.0 

930—1600 

196,4 

H«a,+H«(*)-2Hta 

2.0 

550 — MO 

70.4 

H*Bf,+H«(*)22HtBr 

2.0 

520-950 

63.7 

Sod,  fSn<*)~  2Scd 

2.0 

760-1500 

208,0 

St\Br,+So<M);;2So3f 

2.0 

550-1500 

178,2 

SnJ,  fSn<*)^2Snl 

2.0 

650-1000 

— 

Fe(C0),^Fc+500 

s.o 

350-1200 

68,7 

Ni«O),rNH400 

4.0 

320-1000 

70,0 

Sna,+3Sii<*)-4Saa 

4.0 

BOO— 1200 

425,0 

SnBf,+3Sn(m);t4SnBr 

4.0 

500-1200 

372,0 

Snl4+3Sn(*)^4Sol 

4.0 

500— 1200  • 

— 

HfcOursHea, 

2.0 

— 

1.9 

2.0 

— 

l.t 

2.0 

— 

1.9 

(HF),-6HF 

6.0 

— 

• 

4.0 

700-1200 

—  - 

KEY:  (1)  Dissociating  systems,  (2)  Gas- 

constant  increase  factor,  (3)  Temperature 
of  the  reaction  range,  °K,  (4)  Thermal 
effect  of  chemical  reaction,  kcal/mole. 


The  authors  would  like  to  express  their  thanks  to  mathematicians/ 
programmers  I.  S.  Zakharova,  V.  N.  Pisarchik,  G.  A.  Razumova,  and 
Engineer  Ye.  N.  Bunin  for  the  vast  number  of  calculations  performed. 
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CHAPTER  I 

THE  THERMOOYNAMIC  PROPERTIES  OF 
CHEMICALLY  REACTING  SYSTEMS 


1.  CALCULATING  THE  THERMODYNAMIC 
PROPERTIES  OF  MIXTURES  OF 
REAL  GASES  . 

Let  us  examine  the  following  basic  approaches  to  calculation 
of  the  thermodynamic  properties  of  gas  mixtures  [13].  Here  we  use 
the  term  "partial  pressure"  of  the  gas  in  the  mixture  to  designate 
the  product  of  its  molar  fraction  in  the  mixture  times  the  total 
pressure  of  the  mixture;  we  use  the  term  "equilibrium  pressure"  as 
the  pressure  of  this  gas  in  isothermal  equilibrium  with  a  mixture 
through  a  diaphragm  which  can  be  penetrated  only  by  thiB  gas. 

The  ideal -gas  approach.  In  accordance  with  these  definitions, 
in  the  ideal-gas  approach  we  assume  that  1)  the  equilibrium  pressure 
of  the  gas  in  the  mixture  is  equal  to  its  partial  oressure,  and  2) 
the  equation  of  state  of  the  mixture  is  as  follow  pV  ■  En^T.  In 
place  of  the  second  assumption  we  can  simply  require  that  the  internal 
energy  of  the  mixture  component  be  a  function  only  of  temperature. 

The  relationships  for  the  thermodynamic  properties,  obtained  using 
this  approach,  are  well-known. 

Now  let  us  derive  formulas  for  calculating  the  thermodynamic 
properties  of  pure  gases  using  a  generalized  method  of  passage  to 
the  limit,  since  these  formulas  are  used  In  the  other  approaches. 
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When  deriving  the  formulas,  let  us  direct  our  attention  only 
to  enthalpy  and  entropy. 

The  generalized  method  of  passage  to  the  limit  for  pure  gases. 

Let  us  examine  two  regions  of  pressure: 

2 

a)  low  pressure,  in  which  terms  of  the  order  of  p  are  negligibly 
small  compared  with  terms  of  the  order  of  p  or  one; 

b)  very  low  pressures,  in  which  terms  of  the  order  of  p  are 
negligibly  slight  compared  with  one. 

Let  us  assume  that  at  all  temperatures  the  product  of  volume 
times  pressure  for  a  fixed  mass  of  a  real  gas  can  be  represented  by 
the  following  expression  with  constant  temperature  in  the  region  of 
low  pressures: 


pV  =  F(T,  m)  +  A(T,  m)p  +  0(p*).  (1.1) 

where  F(T,  m)  and  A(T,  m)  depend  for  the  given  gas  only  on  tempera¬ 
ture  and  mass,  while  A(T,  m)  and  its  derivatives  with  respect  to 

2  2 

temperature  are  limited;  0(p  )  are  terms  of  the  order  of  p  .  Further 
let  us  assume  that  at  all  temperatures  the  enthalpy  of  a  fixed  mass 
of  real  gas  can  be  represented  by  the  following  expression  with 
constant  temperature  in  the  low-pressure  region: 

H-HT,  m)+B(T,  m)p  +  0(p*).  '  (1.2) 

where  f(T,  m)  and  B(T,  m)  depend  for  the  given  gas  only  on  the 
temperature  and  mass,  while  B  and  its  temperature  derivatives  are 
limited. 

Substituting  equations  (1.1)  and  (1.2)  into  the  relationship 

we  get 
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s+°*>-7[M§).]  + 

+hr(f).]+ow 

This  equation  is  valid  provided 

Integration  of  this  relationship  with  constant  mass  gives 


F  —  r(m)T,  (1.4) 

where  r(m)  depends  only  on  the  mass  of  the  given  gas. 

Considering  this,  we  obtain  the  following  from  equation  (1.1): 

pV  =  r(m)T  +  A(T,  m)p  +  0(p»),  (1.5) 

while  for  the  region  of  very  low  pressures 

p*V*/rT  « 1.  (1.6) 

Since  volume  is  an  extensive  magnitude,  if  we  replace  V  in  equation 
(1.6)  by  V/m,  r  becomes  a  magnitude  characteristic  for  each  gas. 

If,  however,  the  mass  of  the  given  gas  is  expressed  in  moles,  r  will 
not  depend  on  the  nature  of  the  gas  and  we  will  designate  it  by  R. 
Then 


PV  «  nRT  +  nfipiRT  +  0(p*),  (1.7) 

where  n  is  the  number  of  moles  in  volume  V,  0  is  a  pure  temperature 
function. 

When  integrating  the  differential  relationships  for  enthalpy 
and  entropy  let  us  select  the  following  path  from  the  point  of  the 
normal  state  (pQ,  TQ)  to  an  arbitrary  point  (p,  T),  changing: 

1)  the  pressure  from  Pq  to  p*  with  the  temperature  constant  and 
equal  to  Tgj 

2)  the  temperature  from  TQ  to  T  with  constant  pressure  p# ; 
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3)  the  pressure  from  p*  to  p  with  constant  temperature  T. 


Prom  the  ordinary  differential  relationships  we  get 


»'  r 

H(p,  T) — h (p„  r#)=  J  f nc-/tr+ 

S(p,  T)  —  S(pt,  Tt)  -  -  j*  (-^  )  dp  +  J  -^-dT- 

h*  ”  p*t, 

.  •  -  -m* 


(1.8) 


(1.9) 


In  order  that  the  definite  pressure  integrals  in  formula  (1.9)  not 
diverge  as  the  limit  of  p#  tends  toward  zero,  let  us  add  the  term 
nR/p  to  the  integrand  and  subtract  nR  In  p/pQ  from  all  terms.  We  get 


T*, 


rn/? 

<  dV  \ 

1  p 

^Ip. 

Let 


*  i. 

r 

A*  -  hi  +  J  c^T;  -  H  (/>..  Tt)  + 


(1.10) 


(1.11) 


T+. 

T 


s’-sS+j*  “r'dr:  nS«==s(*-7V>+ 


(1.12) 


4 


where 


r 

+  j '  cvdT,  4- ft 
S>  =  S?+| 


—  RT* 


dT  +  R\nplpr 


(1.18) 


The  Glbbs-Dal ton  law.  The  Gibbs-Dalton  law  states  that  the 
total  pressure  of  a  mixture  of  various  gases  Is  equal  to  the  sum  of 
the  pressures  which  these  gases  would  have  if  they  each  existed 
individually,  at  the  mixture  temperature  and  with  those  chemical 
potentials  which  they  have  in  the  mixture.  This  statement  contains 
an  idea  of  equilibrium,  essential  and  necessary  for  general  thermo¬ 
dynamic  studies  of  gas  mixtures.  Obviously,  the  various  gases 
having  those  chemical  potentials  which  they  had  in  the  mixture  and 
being  at  the  mixture  temperature  will  be  in  equilibrium  with  the 
mixture  through  the  diaphragm  which  is  permeable  only  for  the  given 
gas,  while  the  pressure  of  each  such  gas  is  the  equilibrium  pressure. 

Prom  the  Gibbs-Dalton  law  it  follows  that  the  concentration  and 
number  of  moles  of  a  pure  gas  in  equilibrium  with  a  mixture  through 
a  diaphragm  permeable  only  to  it  are  the  same  as  in  the  mixture, 
hence 

(1.19) 

H  -5X.P  (1.20) 

s-SX  (i.2i) 

Equation  (1.1$)  is  the  ordinary  statement  of  the  Dalton  law.  The 
Dalton  law  follows  from  the  Gibbs-Dalton  law,  but  the  inverse  is  not 
true,  since  the  idea  of  equilibrium  is  lacking  in  the  formulation  of 
the  Dalton  law. 
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The  Qibbs-Dalton  law  leads  to  total  solution  of  the  problem  of 
determining  the  thermodynamic  properties  of  mixtures  of  gases  from 
information  cn  the  behavior  of  the  pure  ijases  of  which  the  mixture 
consists . 


Having  substituted  into  equations  (1.20)  and  (1.21)  the 
expressions  for  the  enthalpy  and  entropy  of  pure  gases  (1.17)  and 
(1.18),  we  obtain  the  following  for  the  enthalpy  and  entropy  of  a 
gas  mixture  subjecc  to  the  Olbbs-Dalton  law: 


.lii]  ] 

JL  Vm 

l  ar  IvJ 

(1.22) 


+  nS]n--+n& 


(1.23) 


The  Gibbs-Dalton  law  leads  to  the  following  expression  for  the  law 
of  mass  action: 


tl-M> 

+  VT-7*)]-  U-25) 

Calculation  of  Kp  for  given  p  and  T  of  the  mixture  is  performed  by 
the  method  of  successive  approximations.  First  we  give  Kp  ■  K*  and 
determine  the  composition  of  the  mixture.  The  mixture  composition 
allows  us  to  estimate  the  right-hand  side  of  the  equation  (1.2*0  and 
obtain  a  second  approximation  of  Kp.  The  process  is  repeated  until 
the  composition  of  the  mixture  becomes  constant. 

The  Lewls-Randal  1  rule.  The  Lewis-Randall  rule  states  that  the 
fugaclty  f^  of  a  gas  in  a  mixture  is  equal  to  the  product  of  its 
molar  fraction  in  the  mixture  times  the  fugaclty  f  of  this  pure 

pH 

gas  with  total  pressure  and  temperature  of  the  mixture: 
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mj\ 


U  amjcJt+ 

The  fugacity  of  a  pure  gas  is  defined  by  the  relationship 

|i-/?Tln/-f  ♦(D.  H»n  —  -  I. 

M  p 


(1.26) 


(1.27) 


Let  us  examine  a  pure  gas  in  Isothermal  equilibrium  with  a  mixture 
through  a  diaphragm  which  is  permeable  only  to  this  gas.  The  chemical 
potential  of  a  gas  l.i  a  mixture  is  p^,  while  that  of  a  pure  gas  is 
p  ,  and  p.  •  u  .  According  to  the  Lewis-Randall  rule  p  and,  conse- 
quently,  ^  can  be  calculated  by  substituting  equation  (1.26)  into 
(1.27) : 


Equation  (1.28)  can  be  rewritten  as  follows: 


(1.28) 


+  RThix,. 


(1.29) 


considering  that  RT  lr.  f  +  (T)  ■  p  .  Let  us  differentiate 

P  K  K  pH 

equation  (1.29)  with  respect  to  p  and  T  under  the  following  conditions 
which  result  from  the  Lewis-Randall  rule:  the  temperature  of  each 
pure  gas  is  equal  to  the  temperature  of  the  mixture;  the  pressure  of 
each  pure  gas  is  equal  to  the  pressure  of  the  mixture;  n.  ■  n  . 

1  n 

We  get 


(1.30) 

(1.31) 


Prom  the  general  theory  of  thermodynamics  we  know  that  the  left 
sides  of  equations  (1.30)  and  (1.31)  are  equal,  respectively,  to  the 
partial  molar  volume  and  the  partial  molar  entropy  taken  with  a  minus 
sign  for  the  1-th  gas  in  the  mixture,  while  the  right  sides  are 
equal  to  the  corresponding  molar  values  for  a  pure  gas. 

Thus, 
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Prom  the  Euler  theorem  on  homogeneous  functions 


(1.34) 

(1.35) 


Using  relationships  from  general  thermodynamics,  we  get 

(1.36) 


Using  equations  (1.15)  and  (1.16),  we  finally  obtain  the  following 
expressions  for  the  enthalpy  and  entropy  of  a  gas  mixture  subject  to 
the  Lewis-Randall  rule: 


(1.36a) 


(1.37) 

However,  the  expression  for  the  law  of  mass  action  is  as  follows: 


(1.38) 


The  generalized  method  of  passage  to  the  limit  for  gas  mixtures. 
The  generalized  method  of  passage  to  the  limit  is  based  on  the 
following  assumption  on  the  behavior  of  gas  mixtures:  at  all  tempera¬ 
tures  the  ratio  of  the  equilibrium  pressure  of  each  gas  in  the 
mixture  to  its  partial  pressure  can  be  represented  in  the  low-pressure 
region  by  the  expression 


£=-«l+D((r.  x.,...)p  +  0(p*), 


(1.39) 


where  depends  only  on  the  temperature  and  composition  of  the 

mixture,  and  its  derivatives  with  respect  to.  T,  r^,  n2>  •••  are 
limited.  Let  un  examine  a  gas  mixture  consisting  of  Zn^  moles  and 
"h"  systems.  Let  us  assume  in  all  measurements  that  the  temperature 
of  each  pure  gas  is  equal  to  the  temperature  of  the  mixture  and  the 
volume  of  each  pure  gas  is  equal  to  the  volume  of  the  mixture,  ^  * 


Let  us  also  assume  that  T,  p^,  p2,  •••  arf-  independent  variables 


for  the  gas  mixture,  T  and  p  ax*e  independent  variables  for  each 

H 


pure  gas.  Then  we  rewrite  equation  (1.39)  as  follows: 


Prom  this 


P«“P*iO  +Oj>). 
2P«“P(‘|  +0mP). 


(1 . 40) 

(1.41) 


where  D  ■  Ex.D. . 
m  11 


Differentiating  equation  (1.41)  first  with  respect  to  one  of 
the  p's  with  T  and  the  remaining  p's  constant,  and  then  with  respect 
to  T  with  all  p ' s  constant,  we  obtain  the  following  for  the  low- 
pressure  region: 


(1.42) 


(1.43) 


Prom  equation  (1.42)  we  have 


-y  -  “  ~y~  0  +  2D«P)| 


while  from  equation  (1.43)  we  derive 


2X.-SU  +  2Dj>). 


(1.44) 


(1.45) 


Dividing  expression  (1.40)  by  equation  (1.44)  we  get 


(1.46) 
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Prom  known  thermodynamic  relationships 


2ff«K“tf(H-2Dmp). 

(1.47) 

For  the  region  of  very  low  pressures  we  have 

p«  =  (28«)*  -p*: 

(1.48) 

.  ^“XC/ 

(1.49) 

s*=2-C 

(1.50) 

h*=2w;. 

(1.5D 

W  -  2/,i/?7'* 

(1.52) 

In  accordance  with  equations  (1.8)  and  (1.9)  t'or  the  enthalpy 
and  entropy  of  a  pure  gas  for  En^  rnolus  of  the  gas  mixture  of 
constant  composition 

p* 

H(p,  T)~H  (pa,  T0)  =  J  dp  + 

Tj>, 

+jc;„r+|[v-r(fjj^  d.53) 

sttn-sw-  J[**-(£y«f+ 


(1.54) 


Considering  expressions  (1.51),  (1.15)  and  (1.48)  we  get 

ff(p*.  r.)  =2//„(p*.  rj «2n.X  =2«.C  (1.55) 

From  equations  (1.50),  (1.16)  and  (1.48)  we  have 

S(p*.  7’c)«2s,ksSL- 

— 2n«*<  * ,n  p»  “  X”*5*  — 2n^ ,n  p,Jfr  (1.56) 


Considering  equations  (1.48)  and  (1.49)  we  get 


(1.57) 


Thus,  from  equalltites  (1.53)  and  (1.55)  we  finally  obtain 

■fcia+j^rfS) J*- 

while  from  equations  (1.54)  and  (1.56) 


(1.58) 


r  i 

\*l\  ] 

1  p 

l  dT  L,\ 

-S(p*.  r.) +Sn, /? In p*= lot,. 


(1.59) 


Let  us  substitute  expressions  (1.57),  (1.58),  and  (1.59)  Into 
formulas  (1 . 53 )-(l.  54) ,  replacing  in  this  cp  ~-  the  limit  p*  by  C, 
considering  equation  (1.46),  and  we  finally  obtain  the  following 
equalities  for  the  enchalpy  and  entropy  of  the  gas  mixture: 


Hu,.  r)-J[v_r(^J<fp+s«lh!, 

S(P,  n- JteHSU.*- 


(1.60) 


(1.61) 


The  expression  for  the  law  of  mass  action,  obtained  by  the  generalized 
method  of  passage  to  the  limit,  is  as  follows: 


(1.62) 


In  view  of  the  fact  that  the  right  side  of  this  equation  depends  on 
the  composition,  it  must  be  calculated  by  the  method  of  successive 
approximations . 

In  order  to  calculate  the  thermodynamic  properties  of  gas 
mixtures  having  constant  composition  using  relationships  derived  by 
the  generalized  method  of  passage  to  the  limit,  we  must  know  the 
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equation  of  state  of  the  mixtu/e  in  such  a  form  that  it  explicitly 
includes  the  number  of  moles  of  the  components.  The  equation  of 
state  of  a  gas  mixture  of  constant  composition  is  usually  given  the 
same  form  as  the  equations  of  state  of  pure  gases,  and  each  parameter 
of  the  mixture,  let  us  say  B  ,  is  associated  with  the  corresponding 
parameters  of  the  components  by  the  relationship  Bm  *  Ix^SjB^j.  The 
parameters  of  interaction  of  unlike  molecules  B^  are  associated  with 
parameters  of  the  components  by  specific  combination  rules  which  have 
been  the  object  of  numerous  theoretical  studies. 

2.  CALCULATING  THE  THERMODYNAMIC  PROPERTIES 
OP  REAL  CHEMICALLY  REACTING  SYSTEMS 


The  thermodynamic  properties  of  real  mixtures  of  constant  com¬ 
position  can  also  be  calculated  from  relationships  derived  by  the 
generalized  method  of  passage  to  the  limit  for  pure  gases,  if  we  know 
the  composition  of  the  mixture,  its  volumetric  behavior,  and  data  on 
the  heat  capacity  of  the  mixture  components  at  low  pressure.  For 
chemically  reacting  real  equilibrium  gas  systems  the  question  of 
calculating  the  thermodynamic  properties  from  the  known  volumetric 
behavior  of  the  systems  (or  from  the  known  equation  of  state)  has  not 
been  developed  with  the  same  degree  of  strictness  as  for  pure  gases 
and  gas  mixtures  of  constant  composition.  In  view  of  the  fact  that 
for  chemically  reacting  systems  the  assumptions  of  the  generalized 
method  of  passage  to  the  limit  for  pure  gases  have  not  been  carried 
out,  the  relationships  derived  for  pure  gases  cannot  be  used  when 
calculating  the  thermodynamic  properties  of  chemically  reacting 
systems.  We  can  show  this  in  another  way.  Let  us  differentiate, 
for  example,  expression  (1.16)  to  calculate  the  entropy  with  respect 
to  pressure,  bearing  in  mind  that  n  ■  En^  Is  a  function  of  the 
pressure  in  terms  of  the  chemical  reaction  equilibrium  constant. 


We  get 


(1.63) 


Since 
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^jainp^O,  (1.6*0 

one  of  the  general  differential  relationships  of  thermodynamics 


is  not  satisfied,  which  proves  the  inapplicability  of  formula  (1.16) 
for  calculating  the  entropy  of  chemically  reacting  systems.  All 
other  methods  for  calculating  the  thermodynamic  properties  of  mixtures 
of  real  gases  having  constant  composition  can  also  be  used  when  cal¬ 
culating  the  thermodynamic  properties  of  chemically  reacting  systems 
with  consideration  of  nonideality.  But  it  is  most  convenient  to  use 
in  this  case  a  method  based  on  the  Lewis-Randall  rule,  for  the 
following  reasons. 

To  calculate  the  thermodynamic  properties  of  chemically  reacting 
systems  with  consideration  of  nonideality  in  all  approaches  we  must 
know  the  heat  capacity  of  the  mixture  components  at  low  pressure  and 
their  volume  behavior,  given  either  in  the  form  of  tables  or  in  the 
form  of  the  equation  of  state.  Ordinarily  the  equations  of  state  of 
mixture  components  are  given  either  for  a  very  narrow  range  or  not 
at  all.  In  this  case,  to  describe  the  volume  behavior  of  the  com¬ 
ponents  we  use  the  law  of  corresponding  state.  When  calculating  the 
thermodynamic  properties  of  substances  which  have  been  studied  only 
little  or  not  at  all,  we  often  use  generalized  tables  of  the  depen¬ 
dence  of  the  fugacity  coefficient  on  the  reduced  parameters.  When 
calculating  the  thermodynamic  properties  of  chemically  reacting 
systems  by  the  Lewis-Randall  rule,  corrections  for  the  nonideality 
of  the  mixture  components  are  very  simply  expressed  in  terms  of  the 
fugacity  coefficients.  It  is  also  easy  to  take  into  account  the 
influence  of  the  nonideality  of  the  components  of  the  mixture  on  its 
composition.  In  this  case  the  composition  is  determined  directly, 
not  by  the  method  of  successive  approximations  as  in  the  case  of 
calculations  using  the  Gibbs-Dalton  law  or  by  the  generalized  method 
of  passage  to  the  limit. 


lH 


The  hypothesis  of  corresponding  states,  advanced  by  van  der  Waals, 
has  been  carefully  discussed  by  Hougen  and  Watson  [14],  They  showed 
that  the  maximum  error  in  predicting  the  volume  behavior  of  a  sub¬ 
stance  is  15^,  while  the  error  in  predicting  the  other  properties  is 
in  the  range  5-35% .  This  shows  the  degree  to  which  real  substances 
are  subject  to  the  hypothesis  of  corresponding  states.  Pitzer  [15] 
has  proposed  another  scheme  for  predicting  properties  which  is  some¬ 
what  more  complex,  since  it  includes,  in  addition  to  the  critical 
parameter,  one  more  parameter,  but  it  gives  results  which  are  one 
order  of  magnitude  more  accurate. 

Statistical  theory  has  shown  [16]  that  a  group  of  substances  will 
be  subject  to  the  principle  of  corresponding  states  only  when  their 
intermolecular  potentials  are  identical  in  form,  although  the  depth 
of  the  minimum  of  potential  and  the  minimum  intermolecular  distance 
associated  with  the  critical  temperature  and  the  critical  volume, 
respectively,  are  characteristic  for  each  substance.  The  intra¬ 
molecular  motions  should  be  classical,  i.e.,  quantum  effects  should 
be  disregarded.  The  only  group  of  substances  corresponding  to  these 
requirements  -  heavy  inert  gases  Ar,  Kr,  Xe  -  correspond  precisely 
to  the  law  of  corresponding  states.  These  are  called  simple  gases. 

We  can  expect  that  various  types  of  the  forms  of  molecules  and 
molecular  dipole  moments  will  lead  to  various  deviations  from  the 
macroscopic  properties  of  simple  gases.  It  was  found,  however,  that 
the  relative  second  theoretical  virial  coefficients  of  a  broad  class 
of  the  molecular  type  of  substances  fall  within  a  single  family  of 
curves  which  can  be  characterized  by  one  parameter.  Molecules  with 
a  high  dipole  moment  are  an  exception.  The  substances  making  up  this 
class  are  usually  called  normal  gases.  According  to  the  principle 
of  corresponding  states,  all  substances  should  have  identical  satur¬ 
ation  length  in  relative  coordinates.  Actually,  this  is  not  so. 

These  lines  comprise  one  family  of  straight  lines  which  differ  from 
one  another  in  their  slopes.  It  was  this  slope  which  was  taken  as 
the  characteristic  parameter  in  describing  the  family  of  curves  of 
virial  coefficients  and  as  the  supplementary  parameter  for  a  more 
precise  system  of  correlation  of  the  nonideality  of  behavior  of 


\ 
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substances.  To  this  we  must  add  that  the  saturation  lines  were 
measured  very  accurately. 

The  third  correlation  parameter  u  is  defined  by  the  following 
relationship: 


where  p_/p  is  taken  with  T„/T  „  -  0.7  at  a  point  rather  far  from 

S  K  p  S  K  p 

the  critical  point.  This  formula  was  selected  such  that  w  is  equal 
to  zero  for  the  simple  gases  Ar,  Kr,  Xe  with  simple  spherical  mole¬ 
cules.  Normal  gases  have  small  positive  values  of  tu.  The  parameter 
a)  is  called  the  acentric  factor,  in  view  of  the  fact  that  it  shows 
the  measure  of  the  deviation  of  the  intramolecular  potential  function 
from  the  potential  function  of  simple  spherical  molecules.  The  p-V-t- 
data  of  a  number  of  substances  were  interpolated  by  series  with  equal 
Pr  -  T  and  constructed  as  a  function  of  w.  Such  graphs  show  that  the 
coefficients  of  compressibility  of  substances  are  described  quite 
well  by  a  linear  dependence  of  w:  Z  =  Zw(p„T;)+«xZ<',0>,.7*r).  The  smoothed 
values  of  Z^  and  \  taken  from  these  graphs,  are  given  in  tables 
Z(0)(pr,  Tr)  and  Z^  ^(Pr«  Tr)  in  works  [15,  17,  18],  Prom  these 
values  of  the  coefficient  cf  compressibility  we  calculated  the 
fugacity  coefficients  and  log  log  \  and  also  corrections 

for  nonideality  of  the  enthalpy  of  substances 

corresponding  to  the  dependences 

log  y  =*  log/0*  +  ulog y<'> ,  (1.66) 


When  calculating  the  thermodynamic  properties  of  gas  mixtures  we  can 
use  the  very  useful  Kay  postulate  [19]  which  states  that  a  gas 
mixture  conforms  to  the  behavior  of  a  hypothetical  pure  substance 
with  pseudocritical  constants,  calculated  on  the  basis  of  known 
critical  constants  of  the  components  of  the  mixture  and  the  composi¬ 
tion  of  the  mixture  from  specific  combination  rules.  Pitzer  and 
Hultgren  [20]  have  proposed  a  generalization  of  this  postulate  with 
consideration  of  the  theory  of  the  acentric  factor  of  pure  substances. 


RT„  j  •  l  RT"  ) 
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3.  CALCULATING  THE  THERMODYNAMIC 
PROPERTIES  OP  A  REAL  CHEMICALLY 
REACTING  SYSTEM 

(N.O^NO.ZJNO+O,;  AlaCI,-2AIC!,;  A!,Bf,22AIBf,) 

Let  us  examine  the  method  of  calculating  the  thermodynamic 
properties  of  real  chemically  reacting  systems  using  as  our  example 
the  mixture  N;043i‘2N02»*2NO+Oj.  A  feature  of  the  behavior  of  the 
chemically  reacting  system  Ns04t*2N0j**2N0+0»,  as  others,  is  the  fact 
that  its  components  N20^  and  N02  do  not  exist  in  pure  form.  They 
always  comprise  an  equilibrium  mixture,  in  which  it  is  impossible  to 
distinguish  the  bonding  of  N02  to  N20^  or  the  dissociation  of  N20^ 
into  N02  from  nonideality  in  their  behavior.  Because  of  this  fact, 
it  is  impossible  to  obtain  individual  data  on  the  volume  behavior  of 
these  components  and,  consequently,  it  is  impossible  to  compile  an 
equation  of  state  for  the  components.  Although  the  volume  behavior 
of  this  system  has  been  studied  rather  well,  particularly  near  the 
exageration  line,  the  absence  of  a  strict  method  for  calculating  the 
thermodynamic  properties  of  chemically  reacting  systems  from  a  known 
equation  of  state  does  not  allow  us  to  use  this  information  for 
calculating  the  thermodynamic  properties.  The  only  way  for  consid¬ 
ering  nonideality  when  calculating  the  thermodynamic  properties  of 
such  systems  is  to  use  some  modification  of  the  law  of  corresponding 
states.  Here,  the  need  for  knowing  the  real  composition  forces  us 
to  ascribe  to  the  components  NgO^  and  N02  specific  critical  parameters 
and,  consequently,  various  saturation  lines.  At  the  same  time,  we 
know  that  the  chemically  reacting  system  N20^  X  2N02  behaves  as  a 
single  substance,  and  not  as  a  mixture  of  two  types  of  molecules, 
having  constant  melting  and  boiling  points,  critical  parameters,  and 
a  single  saturation  line. 

If  we  assume,  for  example,  that  in  the  phase  diagram,  in  addition 
to  the  curve  of  phase  equilibrium  of  the  system  N20^  X  2N02  there 
also  exist  curves  of  the  phase  equilibrium  of  pure  components,  it  is 
easy  to  conclude  as  to  the  possibility  of  achieving  total  dissoci¬ 
ation  of  the  system  by  isothermal  compression  (Pig.  1).  Along  line 
1-2  we  should  have  (3x/8p)T  >  0,  while  with  total  association  of  the 


Fig.  1.  Phase  diagram  of  the 
system  X  2NC2  from  the 

point  of  view  of  the  mixture. 

Designation:  sts  -  atm  (abs). 


molecules  of  the  system  with  lsobario  heating  along  line  2-3, 

(3x/3p)T  <  0.1  Here  x  is  the  percentage  of  dissociated  molecules  of 

n2ok. 

This  conclusion  is  at  complete  odds  with  the  known  behavior  of 
the  degree  of  dissociation  of  the  system  X  2N02  as  a  function 

of  temperature  and  pressure.  Let  us  note  that  the  question  discussed 
in  detail  in  [21]  of  the  possible  deflections  of  the  phase  equilibrium 
curve  says  nothing  in  favor  of  the  existence  of  one,  two,  or  three 
curves,  but  is  connected  merely  with  the  relationship  of  the  rates 
of  dissociation  and  vapor  formation. 

The  critical  parameters  of  the  components  N20ij  and  N02  are 
determined  as  follows.  There  are  two  approaches  toward  defining  the 
critical  parameters  of  substances  which  have  been  little  studied. 

The  first  uses  the  rule  of  additivity,  the  second  uses  experi¬ 
mental  dependences.  The  percentages  of  atoms,  molecules,  ar.d  bonds 
used  in  the  first  approach  are  obtained  from  experimental  data  on 
hydrocarbons,  and  therefore  this  approach  is  Justified  for  organic 
compounds.  It  was  noted,  however,  that  the  percentages  of  N2  and  02 
in  organic  compounds  do  not  correspond  to  the  percentage  of  these 
molecules  in  inorganic  compounds.  From  this  it  follows  that  the 
first  approach  cannot  be  used  for  molecules  of  N02  and  NgO^. 


‘The  mutual  position  of  the  curves  of  the  phase  equilibrium  of 
the  system  and  the  pure  components,  shown  in  Fig.  1,  follows  from 
the  experiments  by  J.  W.  Smith  [22]  on  the  shift  of  equilibrium  with 
drying  toward  the  more  volatile  components. 


We  examined  the  following  methods  based  on  the  second  approach. 
1.  Definition  of  the  critical  temperature: 


a)  the  Smith,  Oreenbaum,  and  Rutledge  method  [23] , 
for  defining  T  : 


_k=lL_ 

(Pi/Pt)**-  I 


+  7»  +  8. 


The  equation 


(1.68) 


where  la  the  density  of  the  liquid  with  T  •  T^,  P2  ie  the  density 
of  the  liquid  with  T  ■  Tj.  Prom  the  EotvSs  and  Segden  relationship 
we  find  that  y'  —  y*  -  r»(|  — 0),/W.  where  6  »  T/Twp,  Y q  is  the  specific 
weight  of  the  liquid  with  T  ■  0.  In  order  to  obtain  suitable  values 
of  T  ,  the  following  conditions  must  be  satisfied: 

Mp 

r,-r,>20.  *3/4r^.  (1.69) 


The  mean  calculation  error  ATKp  •  3 


b)  the  graphic  Altenburg  method  [23).  It  is  necessary  to  know 
the  density  of  the  liquid  at  two  temperatures.  The  mean  definition 
error  AT  •1.5*  (the  maximum  does  not  exceed  5*); 

Mp 

c)  the  Rydell  method  [17].  It  Is  necessary  to  know  the  boiling 
point  at  a  pressure  of  760  mm,  the  temperature  at  which  the  satur¬ 
ated  vapor  pressure  is  100  mm,  and  the  density  of  the  liquid  at  a 
specific  temperature.  The  method  guarantees  definition  of  T  with 

H  p 

a  mean  error  of  0.8S; 

d)  the  Watson  method  [2*1].  According  to  Watson,  for  nonpolar 
substances  capable  of  association  and  dissociation,  we  have  the 
simple  empirical  relationship 

TJT.t  -  0.283 (Wy/-w.  (1.70) 

where  y  1b  the  weight  of  the  liquid  at  the  normal  boiling  point; 

8 

M  is  the  molecular  weight;  Te  is  the  temperature  at  which  the  satur¬ 
ated  vapor  of  1  g-mole  of  liquid  is  contained  in  22.4  l. 
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The  following  equation  hati  been  complied  to  define  T  : 

c 

In  T,  *•  9, a  -Jt  -4.2. 

*1 

Uolng  It,  we  can  construct  a  curve  In  coordinates  Tfl  and  A  -  -  T^, 

with  which,  uolng  the  known  boiling  point,  we  can  uniquely  define  T  . 
Relationship  (1.70)  wao  verified  for  36  nonpolar  substances.  Its 
mean  error  waa  2% . 

Determining  the  critical  density  or  critical  volume.  Tho  moot 
accurate  relationship  for  determining  the  critical  density  Is  the 
following  Ooldhammer  equation  [24]: 

4,  -  4  +  -  o.  (i.7D 

This  relationship  was  checked  for  64  normal  and  abnormal  substances, 
and  good  agreement  with  experimental  values  was  obtained. 

Determining  the  critical  pressure.  The  critical  pressure  can 
be  calculated  using  the  relationship 

■.0,201-0,06*.  (1.72) 

where  w  13  th2  acentric  Pltzer  factor. 

The  value  of  w  can  be  obtained  by  the  Rydell  graphic  method  [17] 
If  this  Is  Impossible,  p  can  be  defined  from  the  empirical  depon- 

K  p 

dence  [24] 

DT 

-7*- -3.75.  (1.73) 

The  value  3.75  is  not  the  same  for  all  substances.  It  depends 
greatly,  as  has  been  established,  on  the  position  of  the  blnodal  In 
the  lg  p-l/T  diagram,  and  It  Is  the  greater,  the  steeper  the  blnodal. 


However,  for  substances  for  which  there  are  no  experimental  data 
these  methods  cannot  be  used. 


To  determine  the  critical  parameters  of  such  substances  we  can 
use  a  relative  method,  l.o.,  use  a  reference  substance  for  which  the 
critical  parameters  are  known.  The  molecular-kinetic  theory  Indicates 
a  connection  between  critical  parameters  and  force  constants  a  (the 
collision  diameter)  and  c  (the  depth  of  the  potential  hole). 


Thodos  and  Plunn  obtained  the  following  empirical  connections 
between  PHp,  THp,  VHp  and  c/k,  o  [25]: 


a  -  0,61  sOi£'; 
a  -  0,561k!/; 
e  -  0,l866klte/\ 

from  which 

•i 

f-M  «,k!5*»M 

from  ZKp  we  define  pKp  by  the  relationship 


(1.7*0 

(1.75) 

(1.76) 

(1.77) 


(1.78) 


(1.79) 


The  collision  diameter  o  can  be  determined  relatively  accurately  from 
data  on  the  crystal  lattice. 


The  critical  temperature  of  N201j  can  be  determined  using  the 
Smith  (a),  Altenburg  (b),  and  Watson  (d)  methods.  We  obtained  the 
following  values:  a)  T  ■  *489 •  3  °K,  b)  T  ■  *489.7  °K,  d)  T  ■ 

Kp  Kp  KP 

-  *495.0  °K.  These  values  agree  within  limits  of  error  guaranteed  by 
each  method.  The  region  In  which  the  errors  of  all  methods  overlap 
gives  a  value  of  the  critical  temperature  THp  ■  491  °K.  The  critical 
volume  of  N^Ojj  Is  defined  from  formula  (1.77) »  the  critical  pressure 
Is  defined  from  relationship  (1.79).  The  critical  parameters  of  N02 
are  defined  by  the  relative  method.  As  the  reference  substance  we 
used  NjO^. 
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NO. 


Thus,  we  obtained  the  following  critical  parameters  of  and 


NA 

NO, 


7^.K  p^,  atm  V».  cm  Vmole 

491  51,6  265,6 

293.6  66,6  94.46 


As  we  know,  the  critical  parameters  of  the  mixture  NgOjj  *  2N02  are 
as  follows  [2]:  T  -  431 . 35  °K;  p  ■  99.99  atm  (teoh);  V  - 

*3  KH 

•  163  cnr/molc.  With  calculation  of  the  composition  of  the  mixture 
N,O,**2N0^’2NO t with  consideration  of  nonideality,  the  attributing 
of  vuriouo  critical  parameters  of  the  components  N20^  and  II02  has  the 
result  that,  in  the  temperature  region  43 1  °K  <  T  <  ^91  °K,  near  the 
saturation  line  of  the  mixture  N20i<  Z  2N02  In  its  gas  phase,  we  obtain 
corrections  for  nonideality  for  the  component  N.,0^  by  the  law  of 
corresponding  states  in  the  liquid  phase,  while  for  the  component  N02 
In  the  gas  phase  this  is  far  from  the  generalized  saturation  line. 

Por  this  reason,  great  error  is  introduced  when  calculating  the  com¬ 
position  of  the  mixture  in  this  region.  A  similar  error  In  this 
region  is  introduced,  both  in  the  enthalpy  and  the  entropy  of  the 
mixture  when  calculating  the  corrections  for  nonideality,  into  the 
entropy  and  enthalpy  of  the  components  and  N02  with  respect  to 

their  pseudocritlcal  parameters.  All  this  greatly  distorts  the  H-S 
diagram  In  this  region.  Less  distortion  of  the  H-S  diagram  near  the 
saturation  line  can  be  obtained  if  w«  calculate  the  correction  for 
nonideality  Introduced  into  the  entropy  and  enthalpy  of  the  mixture 
as  follows. 


Por  each  value  of  temperature  and  pressure  let  us  represent  the 
mixture  N^ogNO^tfNO+O,  by  a  homogeneous  substance  with  pseudocritlcal 
parameters  calculated  by  the  Kay  rule,  and  molecular  weight  equal  to 
M  —  where  MQ  is  the  molecular  weight  of  one  mole  of  the  initial 

nondissoclated  composition  of  the  mixture,  while  i3  the  number 
of  moles  of  the  mixture.  Por  such  a  homogeneous  substance  with 
variable  molecular  weight  and  current  pseudocritlcal  parameters,  in 
view  of  the  change  in  composition,  the  correction  in  the  entropy  and 
enthalpy  of  the  mixture  for  nonideality,  sought  from  the  pseudo¬ 
critlcal  parameters,  lieu  only  in  the  gas  phase  for  temperatures 


431  <  T  <  491  °K,  Although  the  above-indicated  errors  are  preserved 
In  the  composition,  as  before. 

Such  a  calculation  scheme  Is  Identical  to  calculation  under  the 
assumption  of  Identical  corrections  for  nonideality  In  the  entropy 
and  enthalpy  of  all  components  of  the  mixture,  provided  that  these 
corrections  are  sought  from  Its  pseudocrltlcal  parameters. 

Below  we  give  an  algorithm  for  calculating  the  thermodynamic 
properties  of  a  chemically  reacting  real  mixture  NlO,»*2l'KV*2NO+0^  from 
the  scheme  of  a  homogeneous  substance  with  variable  molecular  weight 
and  current  pseudocrltlcal  parameters.  The  law  of  corresponding 
states  Is  used  In  the  Pltzer  modification.  Using  an  analogous  method 
we  calculated  the  thermodynamic  properties  of  the  chemically  reacting 
systems  AI,Clrta2Aia,I  AI,Br,**2AIBrr 


The  equation  for  the  equilibrium  constant  from  the  first  stage 
of  the  reaction  Is  taken  from  [26];  for  the  second  reaction  stage  It 
13  taken  from  [28].  Smoothed  Pltzer  tables  of  log  log  yO» 

(*» 

,  and 

of  and  N02  were  previously  determined;  the  acentric  factors  of 

N20Jl  and  N02  were  defined  from  the  condition  of  agreement  with  the 
Pltzer  relationship: 


— — -V '  are  given  In  [17,  18],  The  critical  parameters 
/ 


The  critical  parameters  of  NO  and  02  were  taken  from  [28];  the 

acentric  factors  were  determined  from  data  on  p  -T„  given  in  [28]. 

8  8 


In  accordance  with  the  assumed  dependence  of  the  equilibrium 
constant  of  the  first  stage  of  the  reaction  on  the  temperature  for 
AH^  (reaction  heat)  and  for  Asj^  (change  in  entropy  during  the 

reaction)  with  T  ■  298.15  °K  we  obtain  the  following  values:  Ah2,  ■ 

0  MA 
■  13,680  cal/mole,  ASj^  ■  42.06  cal/mole*deg.  For  the  second  stage 
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we  obtain  the  following  valuec:  AH82  "  2?,2Q0  cal/mole,  Asp2  " 

■  34.87  cal/mole *deg.  For  we  assume  18  ■  2140  cal/mole 

and  S°98  l5  -  72.72  [293  .  From  this,  for  N02  AH°98  15  -  7910  [293, 

S298.15  “  57l39i  for  N0  ah°98.15  "  21»510  cal/mole,  S298l5  -  50.315. 

For  02  S298  ■  49.0065  [313.  From  the  polynomial  for  the  dependence 

of  heat  capacity  of  N20^  on  temperature  T  [8]  we  obtained  the  below- 
given  polynomials  fcr  the  enthalpy  and  entropy  of  the  component 

The  polynomials  for  the  dependence  of  enthalpy  and  entropy  of  N02, 

No  and  02  components  on  temperature  are  as  follows  [313: 

1)  In  Kpu  --»  14,19712  +  1,121226  Inf-  1,813894  I0-»T  + 

-+0.I31 186- 10— •  7^ _ M7727-IO*  6*03*296  . 

T*  T 


2 )  In -0,5845 In(l0-*r)  +  — _ 

+  17,5121- 0,3079. 10~«r-  1,6017-  10~#ri  —0,1 183. 10~ “7^; 


3)  we  are  given  the  tables 

>!-logv<*>;  fl  — logY*'’:  n‘.r«»T,-TlTui  P,~  /»//>.„ 

where  /«|,  2,  3,  4;  l-N,04.  2- NO,.  3- NO,  4-0,; 

l) )  !*+■»*». 

5)  Xrl=-^p-; 

6) 

7)  <4/&+  /C,17CM-p/C,l)a}-24/Cii}+  (BOtfr-M^KJa}  + 

+  (2/C,i/(>t-80C)<4+  60XW-24/Cia,+4ACi-0; 
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«.-H  |/a}+4  [4p(l-«,y»/iC,t+  <«.+!)) 

2[4p(l-aJ‘//C>1+(a,+  I)j 


9 )  lt  -  37,077  (7  —  298)  +  38,938  \(T*(T*  —  296*)  -  88  78,84 
In/7'298— 84 1,424  •  10*  (I /T—  1/298); 

10)  S| ->34,077 In 77298+73,878- I  (7"*  (7  298) +8878^4(1/7  —1/298)— 320,712-  IO*(l/7^ — 1/298*); 

11)  i,~  10»(/i*+  A_,x-«+A,x  +  V+  V*+  hkx*+hS+hS+  A/+  V*).  x  — 10"*7; 


12)  S|->  A|  In  x+  -=1  X~>+2Atr +—— —  6tx*  + 


3  . 


+  -j- A**4  +  V+  Y  V.  <  3.  4  ( Table  2 ) ; 

13)  x, _ 1--ai— , - _ fcft . ;  x._i. 

I +01+0,0,  **  I  +  O.+  OA*  •  1  +  04+0,0,  t+  2  r 

«  4  4 

1^  )  “  ^^jTcli  Pt »*■  *•»*■  ^*l**|i 


15)  7/j“AW?+ 


16)  S,— 5,+  S?; 

17  )  Aff  - 1.98728  7„(C  +  4i^)(l  +  o,+  0,04); 

18)  AS- ^  +4, 57588 (>1  +  4^5) (1  +  a,+a,04)  fro.  7,--^;  /»,-  JL. 

'  *  4»  Pc* 

19)  //.-  •^•((I-C4)//,+2o,(1— o,)W,+2o,o^/,+o,o//,— AffJ+100; 

2°)  5.--l-((l-a,)S,+2o,(l-e4)S,+  2a,oA  +  o,«4S,-ASJ-l^-I(l-a0ln(l-oJ  + 
+  2a,  ( 1 — a,)  In  2a,  ( 1 — 04) +2a,a,+  lna,o,  +  a,o.  In  o,a,—  (1  +0, + a, <4)  ln(l  +  o,  +  a,og/p]; 


21)  C„ 
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Table  2.  Coefficients  for  calcu¬ 
lating  the  properties  of  the  mixture 
mixture  na22N0,;:2N0+0, 


.. 

f~)  - - 

06mm«h 

NiQt 

NO, 

NO 

Ot 

Th 

491,0 

293,6 

f80,t 

164,6 

Pa 

61 ,00  • 

66,61 

64,00 

60,10 

0.3463, 

0,3740 

0,5767 

0,024 

A//f 

2140 

7910 

21810 

0 

s? 

72,2 

•76,0036 

66,16662 

66,6072 

*• 

— 

-0,24662 

-0.240348 

-0,169746 

hm\ 

— 

4,114. I0'« 

-7,935- 10*« 

-2,669. 10-* 

hi 

— 

6,0/64 

4,83266 

6,20537 

h, 

— 

64,9703 

29,0791? 

30.47837 

h. 

— 

-296  017 

-190,2217 

-146,6616 

h, 

— 

766,17 

276,738 

421,1702 

h* 

— 

-919,213 

-334.5202 

-584,6064 

A, 

— 

44,094 

32.0333 

132.5306 

— 

954.35 

306,769 

625.966 

— 

-649.71 

-212,369 

-400,939 

KEY:  (1)  Arbitrary  designation, 
(2)  Values. 


The  results  of  the  calculation  are  given  In  Table  2  of  the  Appendix. 


4.  CALCULATING  THE  DENSITY  OP  CHEMICALLY 
REACTING  SYSTEMS 

The  density  of  chemically  reacting  gas  systems  is  calculated 
from  the  scheme  of  a  homogeneous  substance  with  veriable  molecular 
weight  by  the  law  of  corresponding  states  in  the  Pitzer  modification. 
The  pseudocritical  parameters  of  a  homogeneous  substance  and  its 
acentric  factor  are  calculated  from  the  following  dependences: 

^er'=  ^*1 7, I*  Pet  "  ^  *Al'»  ®er  “ 

Prom  the  reduced  parameters  Tr=*TlT„,  Pr—PIP„  and  u>cT  we  find 
ZCT:  ZCT=Z‘e»+«i)CTZ<'>.  Prom  this, 

'-^T'  (1-81> 

where  n  -  EnA  Is  the  sum  of  the  moles  of  the  mixture  components,  MQ 
is  the  molecular  weight  of  the  initial  nondissociated  single  mole  of 
the  mixture;  Z^  and  are  tabulated  [17-19]. 
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The  results  of  calculation  of  the  de 
reacting  mixture  N,01^2N012 2NO+0,  are  give 


CHAPTER  II 


METHODS  OF  CALCULATING  THERMODYNAMIC 
CYCLES  FOR  CHEMICALLY  REACTING 
SYSTEMS 


1.  THERMODYNAMIC  CYCLES  FOR  CHEMICALLY 
REACTING  WORKING  MEDIA 

In  most  of  the  stationary  power  plants,  thermal  as  well  as 
atomic,  which  have  been  created  or  are  being  developed,  as  the  working 
medium  for  the  turbines  we  use  steam  with  a  condensation  cycle.  In 
modern  steam-turbine  plants  (STP)  [my]  steam  parameters  have  been 
achieved  relative  to  temperatures  of  565-580°C  and  pressure  2*10  atm 
(abs)  with  a  thermodynamic  efficiency  of  the  condensation  cycle  with 
one  industrial  superheating  of  45-471,  which  corresponds  to  a  41-421 
net  efficiency  of  the  plant  [A],  However,  successes  in  creating 
reactor  units  with  powers  of  600-800  MW  allow  us  to  obtain  an  output 
temperature  of  the  gas  from  type  AGR  reactors  of  650°C,  from  HTGR 
reactors  -  750-850°C.  The  use  of  steam  for  these  temperatures,  as 
shown  by  research  carried  out  at  the  Central  Scientific  Research, 
Planning,  and  Design  Boiler  and  Turbine  Institute  ira.  I.  I.  Polzunov, 
the  Moscow  Power-Engineering  Institute,  the  Siberian  Power-Engineering 
Institute  of  the  Siberian  Branch  of  the  Academy  of  Sciences  of  the 
USSR,  and  others,  is  economically  unfeasible  [1-3] •  Steam  has  a 
second  major  disadvantage  when  creating  high-power  turbine  units: 
the  physicochemical  properties  of  steam  are  such  that  for  the  existing 
structural  materials,  per  unit  discharge  of  the  turbine  we  can  obtain 
powers  to  125-135  MW  [4].  This  creates  considerable  difficulties  in 
designing  large  single-shaft  turbine  units  [4]. 
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In  connection  with  these  features  of  steam,  in  recent  years 
intense  studies  have  been  carried  out  for  both  new  schemes  for  thermo¬ 
dynamic  cycles  as  well  as  for  more  effective  working  media  which  will 
make  it  possible  to  obtain  technical-economic  plant  indicators  which 
are  better  than  in  the  steam-turbine  cycle. 

One  of  the  ways  for  solving  this  problem  can  be  the  use  of 
chemically  reacting  substances  as  the  working  medium.  Certain  of 
these  have  physicochemical  properties  which  are  such  that  they  can  be 
used  as  working  media  for  both  gas-turbine  as  well  as  gas-liquid 
(condensation)  cycles.  For  example,  we  have  the  chemically  reacting 
mixtures  NJO43*2N0I **2NO+Oj,  AtjCl***2AlCI*,  AiaBr»**  2AlBri,  and  others. 

The  possible  use  of  chemically  reacting  (dissociating)  systems 
as  working  media  for  thermodynamic  cycles  and  coolants  was  first 
indicated  in  the  wfcrks  by  Lighthill  (USA)  [9]  and  Kovtun,  Naumov,  et 
al.  (USSR)  [10]. 

Let  us  examine  the  disadvantages  of  a  gas-turbine  cycle  using 
inert  gases.  The  low  effective  efficiency  of  gas-turbine  plants 
(GTP)  [ny]  Is  explained  by  the  high  losses  to  internal  Irreversi¬ 
bility  in  the  cycle.  The  degree  of  influence  of  internal  losses  of 
the  effective  efficiency  of  the  cycle  is  characterized  by  the  work 


ratio 


Lm—  L. 

,==-L  * 


(2.1) 


where  L  is  the  expansion  work,  L  is  the  compression  work, 
m  h 


In  inert-gas  GTP's,  the  work  ratio  Is  low:  'vO.3-0. 5-  This  is 
due  to  the  great  influence  of  losses  in  turbines  and  lines  on  the 
effective  efficiency  of  the  cycle 


(2.2) 


ru  is  the  internal  efficiency  of  the  expansion  process,  n  is  the 

1  K 

internal  efficiency  of  the  compression  process,  Qnofl  is  the  heat 


fed  in  the  cycle.  To  increase  the  work  ratio,  as  follows  from 
formula  (2.1),  we  must  increase  the  expansion  work  or  decrease  the 
compression  work.  Prom  the  formula  for  the  expansion  work 

(2-3> 

we  see  which  parameters  influence  its  value.  Actually,  the  expansion 
work  is  proportional  to  the  gas  constant,  and  increases  with  a 
decrease  in  the  specific  heat  ratio. 

The  compression  work  is  also  proportional  to  the  gas  constant  of 
the  working  medium,  but  decreases  with  decreasing  specific  heat  ratio 

D-L.  (2.4) 

K— 1  T|, 

<  V 

Consequently,  if  the  cycle  is  arranged  so  that  the  process  of 
expansion  occurs  in  the  region  of  a  high  gas  constant,  while  the 
process  of  compression  occurs  in  the  region  of  the  low  gas  constant, 
this  effect  should  lead  to  an  increase  in  the  work  ratio  and,  conse¬ 
quently,  the  efficiency  of  the  plant. 

The  effect  of  a  variable  gas  "constant”  was  first  used  in  a 
cycle  with  supercritical  compression  [29],  where  the  effect  of 
"adhesion"  of  molecules  was  used,  which  lead  to  an  increase  in  the 
molecular  weight  of  the  substance,  while  the  gas  constant  was 
inversely  proportional  to  the  molecular  weight  of  the  substance  R  - 
=  848  M.  The  supercritical  QTP  cycle  was  arranged  such  that  the 
process  of  compression  was  carried  out  in  the  supercritical  region, 
close  to  critical  (Fig.  2).  Then  in  the  compressor  a  gas  having  the 
properties  of  an  Imperfect  gas  is  compressed,  and  the  equation  of 
state  of  such  a  gas  has  the  form 

PV-RTZ,  (2.5) 

where  Z  is  the  coefficient  of  compressibility  which  takes  into  account 
the  difference  between  a  real  gas  and  a  perfect  gas  as  a  result  of 
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Pig.  2.  T-S  diagram  of  gas 
cycles  using  nitrogen  tetroxlde: 
a  -  with  supercritical  com¬ 
pression,  b  -  with  compression 
from  the  critical  point,  c  - 
with  subcrltlcal  compression, 
d  -  critical  point,  (PHp  ■ 

■  99  atm  (abs),  THp  -  431°K); 

I  -  P-1  atm  (abs),  II  -  4, 

III  -  99,  IV  -  110,  V  -  200  atm 
(abs) . 

Designation:  HHan/Hr*rpafl  ■ 

■  kcal/kg'deg. 

the  finite  dimensions  of  the  molecule,  the  presence  of  intermolec Jiar 
forces,  and  association  of  molecules  in  the  supercritical  region. 

For  a  gas-turbine  cycle,  as  shown  above,  it  is  important  that  the 
complex  ZR,  which  can  be  called  the  arbitrary  gas  constant,  be 
significantly  less  during  the  compression  process  than  the  gas 
constant  during  the  expansion  process.  However,  a  great  disadvantage 
of  gas  cycles  with  supercritical  compression  is  the  fact  that  the 
lower  pressure  of  the  cycle  should  be  above  the  critical  pressure  of 
the  selected  gas;  for  most  substances,  however,  it  is  rather  high  and, 
as  a  result,  we  cannot  have  a  high  upper  pressure  of  the  cycle  with 
an  optimum  compression  ratio.  In  addition,  the  compressor  operates 
in  a  region  where  the  gas  constant  changes  very  greatly  as  a  function 
of  temperature,  and  therefore  even  a  slight  change  in  gas  temperature 
at  the  compressor  inlet  (2-5°)  leads  to  an  abrupt  change  in  gas  flow 
through  the  compressor.  Plants  with  such  a  cycle,  operating  under 
variable  conditions,  are  especially  sensitive  to  a  change  in  tempera¬ 
ture. 


This  effect  of  a  change  of  gas  constant  in  a  cycle  with  super¬ 
critical  compression  is  most  clearly  manifested  in  cycles  using,  as 
the  OTP  working  medium,  chemically  reacting  systems  in  which  a  change 
in  the  gas  constant  is  observed  with  thermal  and  radiation  dissoci¬ 
ation  of  the  substances. 


During  the  heating  or  cooling  of  such  systems  there  Is  a  change 
In  the  number  of  moles  and,  consequently,  the  molecular  weight  of  the 
mixture  with  absorption  or  release  of  the  heat  of  ohemlcal  reactions. 

Using  such  systems  as  the  working  medium  In  closed-cycle  OTP's, 
the  gas  in  the  Initial  state  with  maximum  molecular  weight,  l.e., 
minimum  gas  constant,  is  compressed  in  the  compressor  and  heated  in 
the  regenerator  and  heater  to  the  maximum  temperature  of  the  cycle. 
Here  there  is  dissociation  of  the  gas  with  absorption  of  heat  to 
chemical  reactions,  and  an  Increase  in  the  number  of  moles  and,  con¬ 
sequently,  a  decrease  in  the  molecular  weight  and  an  increase  in  the 
gas  "constant."  After  expanding  in  the  turbine  the  gas,  cooled  in 
the  regenerator  and  the  cooler,  recombines  with  the  release  of  the 
heat  of  chemical  reactions  and  a  decrease  of  the  number  of  moles  to 
the  initial  state,  and  the  cycle  is  repeated.  The  high  gas  "constant" 
of  the  working  medium  in  the  turbine  compared  with  the  gas  "constant" 
of  the  working  medium  in  the  compressor  makes  it  possible  to  decrease 
the  percentage  of  power  expended  on  compressing  the  gas  in  the  com¬ 
pressor  to  30-40$  Instead  of  70-80$  for  inert  gases,  which  leads  to 
an  increase  in  the  effective  efficiency  of  the  cycle  using  chemically 
reacting  mixtures  compared  with  cycles  using  inert  gases. 

The  nonmonotonic  nature  of  the  change  in  thermophysical  proper¬ 
ties  of  chemically  reaotlng  systems  considerably  complicates  calcu¬ 
lation  of  the  characteristics  of  the  thermodynamic  cycle:  the 
specific  heat  ratio,  the  heat  of  regeneration,  and  the  regenerator 
effectiveness.  Below  we  examine  methods  for  calculating  the  specific 
heat  ratio  of  a  chemically  reacting  system  and  the  regenerator 
effectiveness  of  thermodynamic  cycles  using  chemically  reacting 
systems  as  the  working  medium. 

2.  THE  SPECIPIC  HEAT  RATIO  OP  A  MIXTURE 
OP  CHEMICALLY  REACTING  OASES 

On  this  question  we  have  the  vast  work  by  Irene  Singer-Bredt 
[ 35 J -  Unfortunately,  although  she  gives  a  correct  general  approach 
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Co  the  problem,  there  are  a  number  of  errors  which  cause  Incorrect 
derivation  of  formulas  for  the  true  isentroplc  indices  for  various 
pairs  of  thermodynamic  Independent  variables. 


In  the  section  dealing  with  chemically  variable  isentroplc  flow, 
with  differentiation  of  the  entropy  of  the  mixture  components  with 
respect  to  Independent  variables  the  fact  Is  disregarded  that  the 
entropy  of  the  mixture  components  Is  selected  not  with  the  total 
pressure  of  the  mixture  but  with  partial  pressures  of  the  components. 
And  since  the  partial  pressure  of  the  component  1*  equal  to  the 
product  of  the  total  pressure  of  the  mixture  times  the  molar  per¬ 
centage  of  the  component  In  the  mixture,  the  derivatives  of  the 
entropy  of  the  component  contain  terms,  derivatives  of  the  molar 
percentage  relative  to  independent  variables,  which  were  not  taken 
Into  account  by  the  author  of  [35] .  This  was  what  led  to  Incorrect 
derivation  of  all  general  formulas  for  the  isentroplc  Indices  of 
chemically  reacting  gases,  and  also  to  the  Incorrect  assertion  that 
for  a  chemically  variable  equilibrium  system  we  do  not  have  the 
equalities 


In  this  work  formulas  were  again  derived  for  the  true  Isentroplc 

Indices  of  chemically  reacting  systems,  and  results  wero  given  for 

calculation  of  the  Isentroplc  Index  In  coordinates  T,  P  for  the 

chemically  reacting  system  N^r  2NQ,-  7NO  +  0^.  Por  comparison,  the 

results  are  given  of  the  calculation  of  tho  equilibrium  specific  heat 

ratio  c  /c 

pc  ve 

Let  us  examine  an  Ideal  gas  In  which  we  have  the  equilibrium 
reaction  »,#,*••  L  h  -  + 

The  equation  for  the  reaction  equilibrium  constant  Kf(T)»  rip?- 

*h®**e  *  ~  Taking  the  logarithm  of  this  equation,  we 

7-i 

obtain  the  following  expression: 
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(2.7) 


13 

Por  convenience  In  further  mathematical  trana  format  Iona ,  let  us 
replace  Xj  by  (  -  the  degree  of  ooapleteneaa  of  the  reaction  [36] 


_  _  "?+  *ft 


where  ~  V a* 


Then  equation  (2.7)  le  written  aa  followa: 


(2.8) 


In  K,(T)  -  V »4  ln(«J  |  —  * In («•■»•  *1)  fvtap. 


(2.9) 


Let  us  derive  the  differential  equations  of  laentropy  for  the  examined 
reaction  In  T,  p;  p,  T;  p,  p-coordlnatea . 

I.  Coordinate  T,  p.  The  entropy  of  the  mixture 


M  ' 


(2.10) 


where  M  la  the  maaa  of  one  mole  of  the  initial  compoaition  of  the 
mixture,  correoponding  to  C“  0;  Sj  la  the  entropy  of  the  i-th  com¬ 
ponent  : 


r.  * 


Since  p^  -  x^p,  while  p^  la  the  atandard  pressure  per  1  atm. 


S,-S?+  ysidT-RtoMj,. 


(2.11) 


The  equation  of  laentropy  with  consideration  of  expression  (2.10)  la 
aa  follows: 


c 


Since  n ,  » n ?  |  v({. 

(»).-■  <  5  ).■  (S)-'K). 

Prom  equation  (2.11)  we  have 

\  tip  Jr  P  Jr,  l  tip  Jr 
From  equality  (2.8)  we  get 


(2.12) 


(2.13) 


(2.14) 


(2.15) 


Let  ua  substitute  expressions  (2.13)-(2.15)  Into  equation  (2.12)  and, 
after  certain  transformations,  considering  that  Xj^  ■  n^/n,  we  get 


0. 


But 


S 


R  (Vjfl*  —  >71?)  R  (>•«*—  M*v) 


(2.16) 


«  0. 


(2.17) 


The  expression  ahead  of  (H/3T)p  and  (H/5p)j  is  transformed  such 
that,  considering  the  solution  (2.7): 


•  •t  t 

^  v(5,  ^  V]  v.  (  5?  I  l  dT )  -  R  v,  In  px, 


XS-R\nK,n, 


(2.18) 
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where  AS  is  the  change  of  entropy  of  the  reaction. 

Using  the  familiar  relationship 

.  „  AS  A  H 

0 

where  AH  is  the  normal  effect  of  the  reaction,  equal  to  A //,- 

-//?  (H^  is  the  enthalpy  of  the  i-th  component  of  the  mixture), 

from  equation  (2.18)  we  get  the  following: 

A H 


•A  - 


(2.19) 


Substituting  expression  (2.19)  Into  equation  (2.16)  and  considering 
the  solution  of  (2.17),  we  get 


(2.20) 


Differentiating  equation  (2.19)  with  respect  to  T  and  making  certain 
transformations,  we  obtain  the  following  expression  for  : 

But,  from  the  Van't  Hoff  Isobar  equation 

\  ir  !r  RT> 


Therefore , 


(£)-«*•(  J^-t) 


(2.21) 


Prom  equality  (2.15)  we  find 


(2.22) 


(ft-' <&~t) 


Substituting  expressions  (2.21)  and  (2.22)  Into  equations  (2.20),  we 
get 


+  -jrf - r-'f - *1  f  -0. 

I  r(^-f)  J' 


(2.23) 


Let  us  show  that  the  expression  In  the  first  set  of  brackets  Is  the 
equilibrium  heat  capacity  cpfi.  For  this  let  us  use  the  equality 


Since 


H 


2* 
M  ’ 


then 


Considering  the  solutions  of  (2.13)  and  (2.21),  we  get 

AH 


(2.2k) 


-CS4-?) 


while  since 


j 
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while 


gv,  if, -Air. 

then 


Thus,  the  final  form  of  the  lsentropy  equation  equation  is  as  follows: 


(2.25) 


II.  Coordinates  p,  T.  The  equation  of  lsentropy,  considering 
that  p  ■  1/V,  looks  as  follows: 


Since 


from  which 


then 


Since 


then 


(2.26) 


(2.27) 
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- -  ma 


ifr.  Tt.— 


(2.28) 


I  dp  Jr  Ml*  T  Mp  dp  Jr 


Considering  expression  (2.13)  and  (2.21)  we  get 


For  ( 3V/3T)p  we  obtain  the  following: 


(2.29) 


Considering  relationships  (2.28)  and  (2.2 9),  after  certain  transfor 
mations  we  finally  obtain  the  isentropy  equation: 


(2.30) 


For  the  examined  thermodynamic  system 


(2.31) 


III 


Coordinates  p(  p.  The  equation  of  isentropy  is  as  follows: 


(2.32) 


but 


/  dS 

(dS\  =  ^  / 

{  dp 

Ir  T  \  dp  Jv' 

l  W  ),  T  ( 

(2.33) 
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with  consideration  of  equation  (2.35) »  the  following: 


Figure  3  shows  the  dependence  of  <p  T  and  cpe/cve  on  temperature  and 
pressure  for  the  mixture  N1Ot*±2NOl»*2NO+Ol . 


As  can  be  seen  from  this  graph,  tc  differs  greatly  from  the 

P  >  * 

specific  heat  ratio  cpe/cve-  Therefore,  in  gas-dynamic  calculations 
and  in  calculations  of  thermodynamic  cycles  with  chemically  reacting 
working  media  when  using  methods  based  on  use  of  the  specific  heat 
ratio  we  must  make  use  of  the  averaged  values  of  Kp  T>  The  nonmono¬ 
tonic  nature  of  the  change  in  Kp  T  with  respect  to  temperature  and 
the  strong  stratification  with  respect  to  pressure  lead  to  great 

errors  when  averaging  k  ™  with  respect  to  temperature  and  pressure. 

P>  * 

This  does  not  allow  us  to  use,  during  the  calculation  of  thermodynamic 
cycles  for  chemically  reacting  systems,  methods  of  calculating  gas- 
turbine  cycles  using  inert  gases  [37]-  The  most  suitable  method  for 
calculating  thermodynamic  cycles  with  such  working  media  is  the  use 
of  the  H-S  (enthalpy-entropy)  diagrams  which,  although  it  iB  laborious 
particularly  (as  will  be  shown  below)  when  calculating  regeneration 
cycles,  involves  no  errors  of  a  numerical  nature. 


3.  CALCULATING  THE  HEAT  OP  REGENERATION 
IN  CYCLES  USING  CHEMICALLY  REACTING 
MIXTURES  WITH  VARIABLE  HEAT  CAPACITY 


Calculation  of  the  regenerated  heat  in  cycles  using  chemically 
reacting  systems  is  complicated  because  of  the  dependence  of  the 
heat  capacity  of  such  systems  on  the  temperature  and  pressure  (Fig.  4), 
which  does  not  allow  us  to  use  ordinary  formulas  to  determine  the 
regenerator  effectiveness  and  the  regenerated  heat. 


Pig.  4.  Dependence  of  the  effec¬ 
tive  lsobarlc  heat  capacity  c__  on 

P® 

temperature  and  pressure. 

Designations:  kms/i/ nr  *  rpsA  - 
kcal/kg*deg;  sts  ■  atm  (abs). 

Below  we  present  calculation  of  the  regeneration  heat  in  cycles 
using  working  media  with  variable  heat  capacity  for  the  case  of 
unidirectional  flow  and  backflow.  Here  we  assume  that  the  mass  flows 
of  working  medium  along  the  heated  and  warmed  sides  of  the  regenerator 
are  equal,  and  pressure  losses  can  be  disregarded. 

Unidirectional  flow.  Let  us  examine  the  behavior  of  the  current 
temperature  head 

*~r*-r  (2.39) 

as  a  function  of  the  change  in  temperature  of  the  warm  wall  of  the 
regenerator,  which  we  took  as  an  independent  variable.  For  this  we 
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use  the  connection  between  the  current  values  of  the  temperature  of 
the  warmed  and  heated  walls  of  the  regenerator  in  terms  of  the  heat- 
balance  equation  for  the  current  value  of  the  temperature  of  the 


heated  wall  of  the  regenerator 


rj-HiPv  r*)-//(pr  T)—H  (p,,  r,). 

For  the  derivative  of  the  temperature  head  with  respect  to  T* 

«'r.  «  l  -  I  re¬ 

considering  equation  (2.40)  and  the  known  relationship  c  ■ 

pe 

we  obtain  the  following  expression: 


(2.40) 


\  or 


c„(ps.  T) 


(2.41) 


From  this  we  see  that  in  the  interval  of  change  of  T*[T^p,  T^],  due 
to  the  fact  that  the  heat  capacity  is  a  positive  magnitude,  r1  0. 
Consequently,  there  is  no  extremum  of  6  within  the  interval  [T^,  T^] 
and  it  is  reached  at  the  ends  of  the  interval.  Obviously,  the 


minimum  6  is  reached  with  T* 


Let  us  show  „hat  QpB|_  will  have 


a  maximum  with  <5mln  ■  0.  Let  U3  examine  the  change  in  QpQr  with  a 
change  in  6: 

Qpn  ~  W(pt.  TJ  -  H  [pK,  T*).  (2.42) 

Considering  dependences  (2.41)  and  (2.42),  for  the  derivative  of 
Qper  with  respect  to  6  we  obtain  the  following  expression: 


dQm  dr* 

dT*  di 


Cp(pv  r)c,(pt,  r») 

Op  (p,i  T)  +  Cp  (p,i  T*) 


(2.43) 


Since  [Q  ]I  /  0,  within  the  interval  of  change  of  6  from  its 

P  B I  0 

maximum  value  to  the  minimum,  we  have  the  possible  value  6*0,  with 

no  extremum  of  Qpar-  Consequently,  both  the  maximum  and  minimum  of 

Q  are  achieved  at  the  ends  of  interval  6.  Obviously,  the  maximum 
per 

Q  is  obtained  with  6  ■  0.  Thus,  with  unidirectional  flow  of  a 
per 

working  substance  with  variable  heat  capacity,  as  well  as  a  working 
substance  with  constant  heat  capacity,  the  minimum  temperature  head 
is  achieved  at  the  output  from  the  regenerator,  while  Qper  is  maximum 


at  6mln*  Calculation  of  the  maximum  heat  of  regeneration  Qp#r 
reduces  to  solution  of  the  equation 

tf(p„  r*) -//(*.  r •) -  H(pt,  T •  -imU)-HU>vT,)  ( 2 .  kk ) 

with  6mln  ■  0,  from  which  we  find  T#  and  substitute  it  Into  equation 

(2.112). 

With  given  regenerator  effectiveness  u  of  the  heat  of  regeneration 

(2.H  5) 

We  must  note  that  It  la  not  very  convenient  to  give  the  regenerator 
effectiveness  when  calculating  the  heat  of  regeneration  In  cycles 
using  working  media  with  variable  heat  capacity ,  since  It  Is  necessary 
each  time  to  calculate  the  minimum  temperature  head  In  order  to 
arrive  at  the  technically  realizable  head.  It  is  much  more  convenient 
to  give  immediately  the  minimum  temperature  head  6  ln,  find  T*  from 
equation  ( 2 . 14 M ,  and,  substituting  T*  Into  equation  (2. ^2),  determine 
Ci  .  This  Is  even  more  true  of  calculations  of  the  heat  of  regen- 

P  01 

eration  in  the  case  of  backflow,  which  we  will  now  discuss. 

Backflow.  The  heat  balance  equation  for  the  given  case  Is 
written  as  follows: 

'/(/V  r*)-ff(P ..  T»)r  H(Pt,  Ttp)  —  H(pv  T).  (2.H6) 


Let  us  examine  the  behavior  of  the  current  temperature  head  with 
a  change  of  T*  from  T^p  to  ,  considering  T2p  as  fixed.  For  we 
obtain,  with  consideration  of  equation  ( 2 . 46 ) ,  the  following  expression: 


f,(Pi.  T*) 
c,(Pv  T) 


(2.47) 


From  this  we  see  that  6  can  have  an  extremum  not  only  at  the  ends  of 
the  interval  [T^p,  T^J,  as  in  the  case  of  unidirectional  flow,  but 
also  within  it.  The  points  which  we  suspect  to  be  on  the  extremum  6 
are  determined  from  the  equation 


e,<Pi.  7*»)-c,(p„  7% 
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(2.148) 


whore  T  Is  defined  from  expression  (2.46). 


Let  us  assume  that  In  the  Interval  of  ohange  of  Tc  equation 
(2.48)  has  one  root  T$  and  it  corresponds  to  the  minimum  d.  This 
minimum  d  was  obtained  for  a  fixed  value  of  T2p.  If  we  are  given 
not  T?p  but  dmJn,  from  the  equation 

c9(pv  f)-c,(*.  (2.49) 


we  can  find  the  root  Tg  which  corresponds  to 
can  find  T0p  from  the  equation 


4mln*  and  from  T8  w* 


H{pv  rj— |r,p-«uv  rt-d»j  (2.50) 


Knowing 
‘min'  " 


T2p»  w®  can  calculate  QpBf. 
follows: 


which  corresponds  to  the  given 


Q^HfjXr 


(2.5D 


Let  us  show  that  maximum  Q  „  Is  obtained  If  we  set  d  ,  •  0.  Por 

per  min 

this  let  us  examine  the  change  of  T2p  with  a  change  In  $Bln.  Since 


I4atolrw 


(2.52) 


with  consideration  of  equations  (2.47),  (2.49),  and  (2.50)  we  get 


from  which 


i4j;w 


f*ky  , 

*,(/V  TJ 


,r»^H 


Cr^r  rt?. . 


(2.53) 

(2.54) 


Prom  equation  (2.54)  It  follows  that  T2p  reaches  an  extremum  at 
the  ends  of  the  Interval  of  change  of  ^mjn»  In  this  case  T2p  Is 
maximum  with  $mln  ■  0.  Prom  this,  Qpg(_  Is  maximum  with  dmln  •  0,  as 
seen  from  equation  (2.51)i  since  enthalpy  Is  an  Increasing  function 
of  temperature . 

If  with  given  $mln  the  root  of  equation  (2.49)  TJJ  corresponds 
not  to  the  minimum  d,  but  to  the  maximum,  l.e.,  at  the  ends  of  the 


Interval  (T^p,  )  we  have  A  which  la  leae  than  the  given  ^min,  4^^ 

•uat  tv  given  either  at  the  right  or  the  left  end  of  the  Interval  of 

change  T#,  depending  on  where  6  in  leaat .  It  la  eaay  to  nee  that  the 

given  #aln  »uat  be  aolected  at  that  end  of  the  Interval  (T^  ,  T*) 

for  which  the  heat  of  regeneration  ia  the  alnlaua  of  th*  two  value# 

calculated  froa  the  enda  of  the  Interval.  Let  u«  aaauae  that  A  ,  la 

min 

given  with  T*  •  Tap*  Then 

t^-iua,  tj.  (2.55) 

where  T„p  •  -  4^^.  Aa  we  aee  froa  equation  (2.55),  Qpar  la 

maximum  with  A  .  •  0.  The  aaae  in  true  if  A  ,  la  glvon  with  T#  • 

oln  aln 

*  Tap*  I»»  l  hi  a  cane, 

Qm  ~  TJ  -  W(*.  T ^  (  2 . 1,6  ) 

7*,,  •  T?  *  w 

Equation  (2.49)  with  given  4a^n  eight  have  no  route  within  the 

Interval  (T^,  J .  In  thin  caae  tho  given  4aln  aunt  bo  aelectcd 

only  froa  the  endn  of  the  Interval  {Tftp,  T^].  If  in  the  entire 

Interval  c,(jv  T*~  f**  4aln  aunt  be  given  with  Tf  •  T^p  and 

Qp(tr  la  calculated  froa  formula  (2.56),  while  if  in  the  entire 

Interval  \T,r  r*- m.  ARln  auat  be  given  with  T#  • 

and  Qp#r  la  calculated  froa  froaula  (2.55).  Here  the  maximum  Qp#r 

lo  obtained  in  the  aaae  aanner  ae  when  A  ,  *0. 

aln 

If  equation  (2.49)  in  the  Interval  of  change  of  T*  hue  several 
roota,  and  wo  are  given  one  of  theae  roots  Anin,  froo  the  value  of 
thin  root  and  the  given  4nin  we  determine  T^p,  while  froo  we 
determine  A  for  all  other  rooto.  Comparing  the  obtained  valuee  of 
A  and  4niij  to  one  another,  we  find  the  aaalleot  of  the  A  and  the 
root  of  equation  (2.4 9)  correapondlng  to  thle  A.  With  thle  value 
of  the  root  we  aelect  the  given  ®Bjn*  and  from  expressions  (2. 50) 
and  (2.51)  we  find  Qpor  • 

Thua,  in  the  caae  of  backflow  of  working  media  with  variable 
heat  capacity,  the  minimum  temperature  head  can  alno  be  within  the 


Interval  of  ciiange  of  the  Independent  variable,  unlike  backflow  using 

working  media  with  constant  heat  capacity,  but  the  heat  of  regeneration 

le  also  maximum  with  A  .  •  0.  The  expreaolon  for  the  regenerator 

min  q 

effect lveneen  in  the  general  case  la  written  as  followo:  u  • 

par  per 

where  Qpwr  le  the  heat  of  regeneration  with  given  A^^,  Qp#r  le  the 

heat  of  regeneration  with  A  ,  •  0. 

aln 

There  le  a  graphic  method  for  determining  the  heat  of  regeneration 
In  cycled  with  working  media  with  variable  heat  capacity  uelng  the 
M-T  diagram  of  the  working  medium.  Figure  5  ohowa  an  example  of  the 
low-  and  hlgh-preeeure  Isobars  corresponding  to  the  hot  and  cold 
walls  of  the  regenerator.  Heat  can  be  regenerated  on  Isobar  sections 
ft' -2  and  2'-k  which  are  cut  off  by  isotherms  T?  •  const  and  T^  •  conat 
Moving  within  this  temperature  interval  the  low-pressure  isobar  such 
that  at  one  point  it  touches  the  high-pressure  isobar  (which  Indi¬ 
cates  equality  of  the  temperature  head  to  zero  at  this  point  of  the 
regenerator),  while  at  all  other  points  the  condition  6  >  0  would  be 
observed,  we  obtain  the  suixlaum  quantity  of  heat  which  can  be  thermo¬ 
dynamically  fed  to  the  regenerator,  Qj,p_2p  "  Hbp  “  H2p‘  In  thl® 
case  the  re generator  effectiveness  is  equal  to  one.  The  regenerated 
heat,  corresponding  to  the  given  regenerator  efficiency,  Q  • 

pur 

•  low-pressure  Isobar  should  be  moved  downward 

(within  limits  of  temperatures  Tj  and  Tj,)  such  that  the  regenerated 
heat  Is  equal  to  Q  r ,  after  which  we  can  determine  here  the  tempera¬ 
tures  of  the  gas  at  the  inlet  and  outlet  of  the  regenerator  along  the 
high-  and  low-pressure  sides,  and  the  minimum  temperature  head.  If  we 
are  given  immediately  the  minimum  temperature  head,  we  can  shift  the 
low-pressure  isobar  downward  such  that  we  obtain  the  given  minimum 
temperature  head,  and  using  the  previously  described  procedure  find 
Qpor  and  the  outlet  temperatures  from  the  regenerator. 

There  Is  a  second  method  for  graphically  determining  the  heat 
of  regeneration  using  H-T  or  H-S  diagrams  and  the  graph  of  the 
dependence  of  the  heat  capacity  of  the  working  medium  on  the  tempera¬ 
ture  and  pressure.  Moving  along  the  low-pressure  isobar  on  the 
graphs  cp  •  f(p,  T)  (see  Pig.  *0,  we  reach  a  position  such  that 


I  1 

Lp 


Pig.  5.  H-T  diagram  In  a  regen¬ 
erator  In  cycles  using  working 
media  with  variable  heat 
capacity. 


between  the  low-  and  high-pressure  Isobars  the  distance  along  the 
temperature  uxls  1b  exactly  equal  to  the  given  $mln-  The  temperature 
on  the  low-prosaure  Isobar,  corresponding  to  this  position,  will  be 
(us  we  can  euslly  Bee)  the  root  of  equation  (2.49),  Calculation  of 
Q  and  the  temperatures  at  the  regenerator  outlet  with  the  known 
root  T*  can  be  done  using  H-T  or  H-S  diagrams,  using  equations 
(2.50)  and  (2. 51). 

As  an  example,  let  us  examine  the  calculation  of  the  heat  of 
regeneration  In  a  cycle  using  a  chemically  reacting  gas  system 
NpOjj  ;  2N02  T  2NO  +  02. 

Equation  (2.49),  depending  on  the  magnitude  and  position  of  the 
range  ['i'2,  TJ(  ]  can  have  In  It  one-three  roots  or  It  may  not  have  them. 
Depending  on  this,  calculation  of  the  heat  of  regeneration  Is  divided 
Into  four  cases. 

I.  In  the  range  [T^,  T^]  with  given  6min  equation  (2.49)  has 
three  roots:  T?0’  T*0’  T5o* 

Let  us  assume  that  $mln  corresponds  to  temperature  T|q,  and  find 
the  temperature  heads  which  correspond  to  the  two  other  points  which 
we  suspect  of  being  the  extremum  -  Tjjg  and  Tjjg.  This  can  be  done  as 
follows.  Based  on  $mln  and  TJo*  fron  tlie  e<*uatlon 
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W(/»..  Tfp)  —  h (p,,  rli  —  &rh)  rj  — //(p,i  r «) 


(2.57) 


we  can  rind  T£p,  and  baaed  on  T from  the  equationa 

tf(P,.  7i#)  -  //CP,.  -  //CP..  rt)  HiPv  T*),  (2.58) 

H(pt.  fu)-H(pt.  -lw-l.-// (P..  r.)- //(/..  7w)  (2.59) 

we  can  find  6^  and  4^  which  are.  respectively ,  the  temperature  heada 
at  polnta  TJQ  and  Comparing  4^ ,  4 ^ ,  and  ^mln  among  themaelvea, 

we  find  the  amalleat  of  them  and,  at  a  temperature  correapondlng  to 
the  amalleat  temperature  head,  give  ^m^n* 

Thla  comparlaon  of  temperature  heada  at  polnta  TJq,  T-  ,  and 
T^0  can  be  done  without  aeeklng  ,  4 2,  and  4^.  Let  ua  aubtract 
equation  (2.57)  from  expreaalona  (2.58)  and  (2.59).  We  obtain  the 
following  relatlonahlpa: 

«(P,.  rj»— a,)— //  ov  -&«)-.«  (p,.  r;t)  -  //  (p..  ri). 

r;-a,)-//(p,.  r‘»-a«)  -//(p,.  r:a)-Hu>,.  r»y 

In  these  equalltle8  In  place  of  4?  and  4^  we  ai.batltute  4mlR,  and 
find  the  sign  of  the  following  expression: 

//  (P«.  7k  —  4««)  -  H  (P..  Tut  —  (m)  — ■ 

—  H  (pt.  rid  -  H  (p,,  7*»);  (2.60) 

W  (p,.  r»  —  6r>.»)  —  H  (p,,  T*  —  4«»)  — 

-«(p4.  r*«)  «(p,.  r»).  (2.6D 

If  both  expreaalona  are  greater  than  zero,  obvloualy  the  temperature 
head  at  polnta  Tj|0  and  T|q  la  greater  than  6min  and,  consequently, 
4min  corresponds  to  T|q.  If  at  leuat  one  of  expreaalona  (2.60)  or 
(2.6l)la  leas  than  zero,  the  temperature  head  correapondlng  to  It  la 
leaa  than  ®mln  and,  consequently,  4mln  must  be  given  with  the  value 
of  the  root  of  equation  (2. ^9) .corresponding  to  thla  expression. 


If  both  expreaalona  (2.60)  and  (2.61)  are  leas  than  zero,  both 
42  and  4^  are  leas  than  imln-  In  this  case,  assuming  that  4fflln 
corresponds  to  T|q,  from  the  sign  of  the  expression 


K9 


(2.62) 


//  (P|,  T  m  ~~  bt  «)  ■”  H  IP).  7";t  •) 

-  //ip*.  r!o  -r  #i/v  rk# 

wo  rind  the  temperature  head  at  point  T|q.  If  expression  (2.62)  Is 
greater  than  zero,  the  temperature  head  at  point  TJjg  Is  greater  than 
at  point  Tj50  and,  consequently,  ^mln  oust  be  given  at  point  T|Q.  If 
expression  (2.62)  Is  less  than  zero,  the  temperature  head  at  point 
T^0  Is  less  than  at  point  T|q,  and  6mjn  must  be  given  at  TJJq.  Equal¬ 
ity  of  expressions  (2. 60)-( 2 . 62)  to  zero  obviously  indicates  equality 
of  the  compared  temperature  heads.  Then  6mln  can  be  given  at  any  of 
the  compared  points. 


Having  defined  the  root  of  T|Q,  at  which  we  must  give  6mjn» 
from  the  equation  Hip,.  T.r\—  H(p..f»  —  6m^  mH(ot.Tt)— Hto,.  T\)  we  find  T2p 
and  then  calculate  the  heat  of  regeneration  Qm  -  //(p,,  Ttp)~  H(pv  TJ. 

If  we  seek  the  roots  of  equation  (2.49)  In  ascending  order,  i.e., 
T»<T’»<  T»  as  we  see  from  the  graph  of  the  dependence  of  the  heat 
capacity  of  the  chemically  reacting  system  N/>,«<k2NOl«*2NO  O,  on  the 
temperature  and  pressure  (see  Pig.  4),  by  replacing  the  sign  of 
expression  ct(pv  T*)  —  cp(pt.  T*  —  on  passage  through  these  roots  we 

can  conclude  that  and  will  always  correspond  to  minimum  6, 
while  T*Q  corresponds  to  maximum  6.  Prom  this  It  follows  that  minimum 
6  is  never  reached  at  the  ends  of  the  interval  T#[T2,  T^],  and  we  can 
compare  the  temperature  heads  only  at  the  two  points  corresponding 
to  Tf0  and  T*q. 

II.  In  the  range  [T2,  T^]  equation  (2.49)  has  two  roots  - 
and  T«0. 

Using  the  same  method  as  in  the  case  of  three  roots,  let  us 
compare  the  temperature  heads  at  the  points  corresponding  to  the 
roots.  Let  us  compare  the  smallest  of  these  with  the  temperature 
heads  at  the  ends  of  the  interval  [T2,  T^]. 

In  order  to  compare,  for  example,  the  temperature  head  at  the 
end  of  with  the  temperature  head  at  TJq,  it  Is  necessary  to 
determine  the  sign  of  the  expression  H{pv  —  — H[pv  7*I§  —  An(s)  —  H (/>,, 
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7*«)  -f  H (pt.  T'n),  which  ia  obtained  from  the  same  concepts  as  expressions 
(2.60  M2.62). 

III.  In  the  interval  [T2,  T^]  equation  (2.49)  has  one  root  Tjp. 

This  root  can  correspond  either  to  the  minimum  or  the  maximum 
of  the  temperature  head.  This  can  be  determined  by  alternating  the 
sign  of  upon  passage  of  T*  through  TJQ,  i.e.,  by  changing  the 
sign  of  the  expression 

<V(/V  T»).  (2.63) 

If  TJ0  corresponds  to  the  minimum  temperature  head,  by  setting 
it  equal  to  $mln  we  can  define  Qper*  If  T*Q  corresponds  to  the 
maximum  temperature  head,  the  minimum,  consequently,  is  found  at 
either  end  of  the  Interval  [T2,  T^].  This  is  easy  to  establish  by 
determining  the  sign  of  the  following  expression: 


H{pv  T,  f  8.J- 


-//(/>,.  Tt-6mJ  +  H(p„  rj. 


(2.64) 


If  this  expression  is  less  than  zero,  the  temperature  head  is 
obviously  minimum  at  temperature  Tj.  In  this  case  the  heat  of  regen¬ 
eration  is  calculated  as  follows:  (P«.  TJ-H  </>*.  Tt  + 

If  expression  (2.64)  is  greater  than  zero,  the  temperature  head 
is  minimum  at  temperature  and  the  heat  of  regeneration  Qpar* 
h(p2,  -  «mln)  -  h(p2,  t2). 

IV.  Equation  (2.49)  has  no  roots  in  the  interval  [T2,  T^]. 

Prom  this  it  follows  that  the  minimum  and  maximum  temperature 
heads  are  located  at  the  ends  of  the  interval  [T2,  T^].  If  6^,  >  0 
in  the  entire  Interval  [T2,  T^],  i.e., 

c,(Pv  r*)>0,  (2.65) 
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the  minimum  temperature  head  obviously  corresponds  to  T2,  and  the 
heat  of  regeneration  -  «(p4,  r,  +  6,J. 

If  6^  <  0  in  the  entire  interval  [T2,  T^],  i.e., 

ePU>t,  --r,(p4.  r*)<0.  (2.66) 

» 

the  temperature  head  is  minimum  at  and  the  heat  of  regeneration 
Qvtt  ’*  H  (P»*  — ®ibi»)  —  W(p,,  7\). 

Based  on  the  above,  we  can  write  as  follows  the  algorithm  for 
calculating  the  heat  of  regeneration  in  the  cycle  with  N20^  t  2N02  t 
t  2N0  +  02. 

I.  Equation  c,(p4,  T*)^cf{pv  T"  —  has  the  roots  TJQ,  T|Q,  TjjQ  and 
T*0  <  T*0  <  T*0- 

1)  If  H(pv  TM-6m.J  —  T\%— 6„J  —  H{pv  7*)  +  H{pt,  T»): >0.  proceed  to 

No.  2;  <  0,  proceed  to  No.  3; 

2)  if  lt(p,,Tl»—&VllJ  —  mpvTn  —  6mlJ  —  mpvT*)+mpt.T\o):>0,  proceed  to 
No.  <  0,  proceed  to  No.  6; 

3)  if  «(p,.  7jo-<W  —  W(P».  T»— ».|J  —  H{pt,  7*)  +  //(p4,  T'm):>0,  proceed  to 
No.  5;  <  0,  proceed  to  No.  6; 

*0  the  equation  for  finding  Tu  H[pv  Tit)  —  H{pv  7*I»  —  °  Hipt.  Tt)  - 

“  H(Pi*>  Tf0)i 

5)  the  equation  for  finding  T2p:  H(p2>  T2p)  -  H(p2,  T|Q  - 

-  W  ■  H(<V  V  -  T5o)l 

6)  the  equation  for  finding  T2p:  H(p2,  T2p)  -  H(p2,  T»Q  - 

-  W  -  V  T?0)J 

7) 
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II.  The  equation  cp(pt,  f*)- ep[pt,  T*— fl^)  has  the  roots  TJQ,  TgQ. 

1)  If  H(pt,Tlt~imin)-H(pt,  Tw-6min)-H{pv  T»)+H(pt.  fio):  >0,  go  to  No.  2 
<0,  go  to  No.  7; 

2)  if  H(pv  Tt- i**)- Hip,,  T‘,»-6miJ-H(pv  TJ  +  Hfa.  T^):  >0.  go  to  No.  3; 
<0,  go  to  No.  5; 

3)  the  equation  for  finding  T2p:  H(p2,  T2p)  -  H(p2,  T»Q  - 

-  W  ■  H“v  V  -  H<>v  T!o); 

■O  0,,-hip,.  T„)-  Htpr  r,y. 

5)  Tt  —  fimlnl 

6)  rs,)-//(ft,  rj; 

7)  if  H(pvT'it-6Mn)-H(pvTt)-H(pi.  ri)+H(p4>  f, +  «*,„):  >0. 

8)  r0-rt  +  w 

9)  Q„~H(pv  TJ-H(plt  T„); 

10)  the  equation  for  finding  T2p:  H(p2,  T2p)  -  H(p2,  T*Q  - 
6min^  “  H^pi|»  Ti|)  "  h^p4*  T?0}J 

11)  Q„r~H(pt,  Tlf)-H{pv  Tt). 

III.  The  equation  «-p(p,,  T,)  =  cp(p..  T*  —  ^)  has  one  root  T«Q. 

1)  If  cp(Pj.  T*—  ftmin)— cp(pt.  T*)  on  passage  through  TJQ  changes  sign 
from  +  to  go  to  No.  2;  from  -  to  +,  go  to  No.  5; 

2)  if  H(p,.  Tt)-HU>v  rt  +  6mIO)-//Cp,.  rl-6min)  +  //(pJ>  T,):  >0.  go  to  No.  3; 
<0,  go  to  No. 

3 )  <?„  -  H  </v  r4  -  **)  -  //  <pr  r,); 
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*0  Qjxt  ="  rt)  —  H(pt,Tt  +  i, nu); 


5)  the  equation  for  finding  T2p:  H(p2  T2p  )  -  Htpj.  t;0  - 

-  W  ■  H(»V  V  '  Ht»H-  T!o>i 

6)  Qxr  *  H{p..  Tip)  —  H(pv  Tt). 

IV.  The  equation  cp(pv  T*  —  6nun)- <•,(/>«.  r»).  No  roots. 

1)  If  T»  —  bmub-c,{pv  r*):  >0,  go  to  No.  2;  <0,  go  to  3i 

2 )  <?,«,  -  W(pt.  TJ-H(pit  r,+w. 

3 )  Qp*.  -  H  (p„  r4  -  6*,.) —H(pt,  TJ. 

For  systems  of  the  type  AljCl^tfAlCI,,  AI,Br^±AIBr,  the  heat  of  regen¬ 
eration  can  be  calculated  much  more  simply  in  view  of  the  fact  that 
equation  (2.^9)  has  only  one  root  or  none.  In  this  case,  if  there 
is  one  root,  it  must  correspond  to  the  minimum  temperature  head,  and 
there  is  no  need  to  verify  the  minimum  of  the  end  of  the  interval 
[T2,  Tj,]. 

Algorithm  for  calculating  the  heat  of  regeneration. 

I.  Equation  cp{pv  T*  —  «min)  =  cp[pv  T*)  has  the  root  Tg: 

1)  the  equation  for  finding  T2p:  H(p2,  T2p)  -  H(p2,  Tg  - 

-  W  ’  H(p4>  V  "  H(IV  T0)! 

2)  Q„=H(p„T„)-mt„TJl. 

II.  The  equation  cp (p,,  T*  —  Asm)  =  cp (pt,  T*)  has  no  roots: 

1)  if  cp{pv  T*— <U.)-Cp(pt,  T*)\  a)  >0,  go  to  No.  2;  b)  <0,  go  to  No.  3; 

2)  Tt)-H(pt.  T.  +  W. 

3)  <?„,,=*  #(/>*,  — ftnia)  — W(P,,  7**). 
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4.  FEATURES  OF  THE  THERMAL  CALCULATION 
OF  HEAT  EXCHANGERS  USING  CHEMICALLY 
REACTING  MIXTURES 

The  calculation  of  heat  circuits  for  power  plants  using  dissoci¬ 
ating  gases  has  certain  features  and  difficulties.  It  is  quite  com¬ 
plex  to  calculate  the  heat  of  regeneration  in  cycles,  the  thermal 
calculations  of  intermediate  heat  exchangers  in  binary  plants,  coolants, 
etc.  These  difficulties  are  caused  by  the  nonlinear  nature  of  the 
change  in  thermodynamic  properties  of  the  dissociating  gases  as  a 
function  of  temperature  and  pressure  [8];  this  in  turn  leads  to  a  non¬ 
montonic  change  in  temperature  head  along  the  heat  exchanger,  as  a 
result  of  which  the  minimum  temperature  head  can  be  found  not  only  at 
the  ends  but  also  in  the  inner  regions  of  the  heat  exchanger.  There¬ 
fore,  thermal  calculation  of  heat  exchangers  must  be  done  with  this 
condition  taken  into  account.  Below  we  give  a  solution  of  the  problem 
of  finding  the  pcint  of  minimum  temperature  head  and  a  method  of 
thermal  calculation  of  a  heat  exchanger  with  a  given  minimum  tempera¬ 
ture  head. 


The  problem  was  solved  for  a  backflow  circuit  of  motion  of 
coolant  without  considering  losses  of  pressure  and  with  different 
flows  along  the  heating  and  warmed  sides  of  the  heat  exchanger. 


Figure  5a  shows  a  diagram  of 
as  a  function  of  heat  flow. 


the  distribution  of  temperatures 


Fig.  5a.  Diagram  of  the  dis¬ 
tribution  of  temperatures  in  a 
heat  exchanger  as  a  function  of 
ht. at  flow. 


The  form  of  solution  of  the  problem  will  depend  on  the  manner  in 
which  the  initial  data  are  given.  We  will  examine  two  methods  here. 
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I.  Let  us  consider  that  the  following  values  are  given  (see 
Pig.  5a):  Tjj,  ,  p2,  px,  Q2,  Qx,  6ml  .  We  must  determine  the 
output  temperatures  of  the  coolant,  l.e.,  T2  and  T^.  We  are  given 
the  value  T2<  Then  the  heat-balance  equation  for  the  i-th  section  Is 
written  as  follows: 

T, -6J-H 0>„rjj.  (2.67) 

where  H(p,  T),  the  enthalpy  of  the  gas.  Is  a  function  of  temperature 
and  pressure;  Q.^  is  the  heat  flow  for  the  i-th  section  of  heat 
exchangers;  Ti  is  the  current  temperature  in  the  i-th  section  along 
the  hot  side  of  the  heat  exchanger;  TA  -  6^^  is  the  current  temperature 
in  the  i-th  section  along  the  cold  wall  of  the  heat  exchanger;  is 
the  current  temperature  head  for  the  i-th  section;  G2  is  the  flow 
of  coolant  along  the  hot  side  of  the  heat  exchanger;  Gx  is  the  flow 
of  coolant  along  the  cold  side  of  the  heat  exchanger. 


Let  us  differentiate  expression  (2.67)  with  respect  to  current 
temperature: 


Since 


while 


(2.68) 


where  cp(p,  T),  the  isobaric  heat  content  of  the  gas,  is  a  function 
of  temperature  and  pressure,  by  substituting  these  transformations 
into  equation  (2.68)  we  get 


T-*) 


(2.69) 
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After  certain  transformations  of  expression  (2.69)  we  obtain 
the  following: 


dfl  =  Q|  Cp  (Pt<  T)  _ I 

VT  0,  cp(px,  T-6), 


(2.70) 


Let  us  verify  function  <5(T)  for  the  extremum,  equating  its 
derivative  to  zero: 


0j_  cB(p„  r»)  =  , 

Gx  Cp(px,  T»-S*)  ’ 


(2.71) 


i.e.,  Gtcp(pt,  T*)  =  Gxcp(px,  T‘  —  6*).  where  T#  is  the  extremal  point;  6*  is 
the  temperature  head  at  the  extremal  point,  equal  in  magnitude  to 

Smln 

As  can  be  seen  from  expression  (2.71),  the  extremal  point  in  a 
heat  exchanger  can  be  reached  when  the  "water”  equivalent  Qcp(p,  T) 
are  equal. 


Thus,  given  the  value  6*  =  6mln,  we  can  always  find  the  extremal 
point  T*  (if  it  exists  in  the  investigated  range  of  temperature). 

The  unknown  values  of  the  output  temperatures  of  the  flow  (T2  and  T^) 
can  be  found  from  the  system  of  equations 


Gx[H(px.  T*  —  b*)—H (px,  7\)] 
=  G,[W(p,.  T*)  —  H(pt,  7V)]. 
Gx[H(Px,  Tt)—I1(px,  T*-6»)] 
=  Gt[H(p,,Tt)-Hipt.  7*)]. 


(2.72) 


II.  We  are  given  the  following  values:  T1,  Tg,  T^,  Gg,  5mln* 

We  must  determine  T,  and  G  . 

3  X 

In  this  case,  solving  the  system  of  equations  written  for  section 


C,  H(px,  7**  — 6*) — H(p„  Tx) 
Gx  H{p..  T')-H(pv  TJ  ' 

G,  c,{px ,  r»-6») 

Gx  c,(pv  T»)  ’ 


(2.73) 


we  obtain  the  following  expression  for  determining  the  position  of 
the  extremal  point  T*: 


H(pt,  T*)—H{p.,  tj  //(p^r»-a«)-//(p^r4> 

cpU>vT*)  ”  c,(p„  r*-a») 


(2.7*0 


Equation  (2.7*0  allows  us  to  find  the  extremal  point  T*  with  given 
temperature  head  in  it,  and  then  determine  the  unknown  flow  and  the 
unknown  output  temperature,  solving  the  following  system  of  equations: 


o  G  *0V  r*)-//(p,,  Tj> 


W(Pr  ?*)  =  H  (Pr  TJ  +  Gt 


H(p„  Tt)  —  H(pt,  7V 


(2.75) 


We  examined  two  basic  methods  for  solving  this  problem.  The 
point  of  minimum  temperature  head  can  be  found  either  from  equation 
(2.71)  or  from  (2.7*0.  Let  us  study  the  solution  of  these  equations. 

1)  Equations  (2.71)  and  (2.74)  can  contain  several  roots,  i.e., 
several  extremal  points  T*.  In  this  case  we  select  that  solution 
which  satisfies  the  condition  6  >  6min  in  the  investigated  temperature 
range . 

2)  The  presence  of  an  extremal  point  in  the  investigated 
temperature  range  assumes  a  minimum  In  it  as  well  r>3  a  maximum 
temperature  head,  i.e.,  in  the  latter  case  6  <  ^m^n»  which  leads  to 
incorrect  calculation  of  the  heat  exchanger.  To  avoid  possible 
errors  in  the  calculation,  we  must  verify  the  boundary  points  at  the 
minimum. 


3)  If  the  flows  of  coolant  are  identical  along  the  heating  and 
heated  walls  of  the  heated  exchanger,  equation  (2.71)  assumes  the 
form  obtained  previously: 

c,(p»  r*)  =  cp[px.  r*  - 6*).  (2.76) 
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CHAPTER  III 


GAS  AND  GAS-DYNAMIC  CYCLES  USING 
CHEMICALLY  REACTING  WORKING  MEDIA 


1.  THE  QAS  CYCLE 

Unlike  gas-turbine  plants  using  inert  gases  (He,  N2>  and  others), 
in  cycles  using  dissociating  gases  we  should  expect  an  influence  of 
the  shape  of  the  boundary  curve  on  the  thermodynamic  indices  of  the 
gas-turbine  cycle,  since  most  of  the  dissociating  systems  pertain  to 
substances  with  relatively  high  boiling  point.  In  addition,  at 
atmospheric  pressure  in  gas  cycles  with  dissociating  gases  we  should 
expect  a  significant  influence  of  the  absolute  pressure  of  the  gas 
on  the  thermodynamic  indices  due  to  the  dependence  of  the  thermo- 
physical  properties  not  only  on  temperature  but  also  on  pressure. 

During  thermodynamic  studies  of  gas  cycles  we  examined  a  number 
of  closed-cycle  GTP's  using  dissociating  N20i<  gas.  We  studied  a 
simple  circuit  (Fig.  6c),  a  circuit  with  regeneration  of  heat 
(Fig.  6a),  anr'  a  circuit  with  regeneration  and  intermediate  heating 
of  the  gas  (Fig.  6b) . 

Thermodynamic  calculations  of  simple  gas  cycles  for  GTP's  using 
dissociating  gas  with  a  change  of  minimum  temperature  of  the 

cycle  in  the  range  iiO-158°C  and  maximum  pressure  of  6-200  atm  (abs) 
and  temperature  of  600-1200°C  show  .d  that  the  optimum  compression 
ratio  7T  is  60-80,  which  considerably  exceeds  the  optimum  values  of 
the  compression  ratio  using  inert  gases  [37]. 
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Pig.  6.  Circuits  for  gas  cycles: 
a  -  with  regeneration,  b  -  with 
regeneration  and  intermediate 
superheating,  c  -  simple. 

KEY:  (1)  gas  generator,  (2)  com¬ 

pressor,  (3)  regeneration,  (4) 
reactor,  (5)  turbine,  (6)  high- 
pressure  turbine,  (7)  low- 
pressure  turbine,  (8)  inter¬ 
mediate  superheater,  (9)  coolant. 


The  occurrence  of  chemical  reactions  with  absorption  and  release 
of  heat  during  heating  and  cooling  has  a  favorable  influence  on  the 
effectiveness  of  regeneration  in  the  cycles,  which  results  in 
increased  cycle  efficiency.  Therefore,  chemically  reacting  working 
media  are  used  to  best  advantage  in  cycles  regeneration. 

When  analyzing  gas  cycles  with  regeneration  we  examined  the 
influence  of  the  minimum  parameters  (temperature  T1  and  pressure  p-^), 
the  maximum  parameters  (temperature  and  pressure  p^),  the  coeffi¬ 
cient  of  low-space  and  high-pressure  losses,  and  the  heat  of  regen¬ 
eration  on  the  effectiveness  of  the  cycle. 
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2.  THE  INPLUENCK  OP  THE  PARAMETERS  OP 
THE  OAS  CYCLE  OH  ITS  EPPICIENCY 


Tht  Influence  of  heet  regeneretlon  on  the  effectiveness  of  the 
cycle.  To  raise  the  thermal  economy  of  OTP's  we  uoe  regeneration  of 
heat  In  the  cycle;  In  this  case  the  heat  of  the  gao  spent  In  the 
turbine  Is  used  to  heat  a  compressed  gas  which  makes  it  possible  to 
decrease  the  irreversibility  of  the  cycle.  An  analysis  of  gas  cycles 
using  N^Ojj  as  the  working  medium  showed  that  the  use  of  regeneration 
increases  the  cycle  efficiency  by  a  factor  of  1.5-2  and,  in  addition, 
shift*  the  optimum  efficiency  toward  lower  valuea  of  the  compression 
ratio  (ti  )  (Fig.  7),  which  decreases  the  dimensions  of  the  compressor 
and  l.e  turbine.  The  shift  of  maximum  efficiency  toward  lower  values 
of  n  is  explained  by  the  fact  that  at  low  values  of  n  the  tempera- 

K  K 

ture  at  the  end  of  lsentroplc  compression  rises  insignificantly; 
thus,  the  region  of  the  Isobar  in  which  regeneration  is  possible 
increases,  approaching  and  T^,  respectively. 


Fig.  7.  Dependence  of  thermal 
efficiency  on  the  compression 
ratio  ith:  T3  -  1073°K;  Tx  * 

*  390°K;  p2  *  50  atm  (abs); 

Va  ’  °'91i  °H.a  ‘  0 * 90 * 
nT  -  0.90;  nK  -  0.88. 


When  examining  regeneration  in  a  cycle  with  nondissociating 
(inert)  gases  it  is  theoretically  possible  to  permit  total  regener¬ 
ation  of  the  heat  in  a  given  temperature  range,  i.e.,  regeneration 
can  be  considered  a  reversible  process  in  the  presence  of  an  infinite 
heat-exchange  surface. 

However,  regeneration  in  a  gas  cycle  with  dissociating  gases 
even  theoretically  is  an  irreversible  process.  Here  the  completeness 
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oT  regeneration  la  eatlmated  In  ratio  to  tho  maximum  posnible  quantity 
of  heat  which  can  be  used  In  the  regeneration,  l.e.,  when  the  minimum 
temperature  head  Is  equal  to  zero.  The  magnitude  of  regenerator 
effectiveness  depends  on  the  value  of  the  minimum  temperature  head  In 
the  regenerator  and  the  temperature  range,  and  In  practice  varies 
from  zero  to  one. 

In  the  case  at  hand  we  examined  the  Influence  of  the  value  of 
minimum  temperature  head  on  the  effectiveness  of  the  cycle.  Figures 
8  and  9  show  the  dependences  ( ^m^n)  and  the  dependence  Qper  <w- 

As  can  be  seen  from  the  graphs,  regeneration  raises  the  Internal 
efficiency  very  abruptly,  with  the  optimum  n4(*H)  shifted  toward  high 
values  of  ir  with  an  increase  in  minimum  temperature  head  in  the 
regenerator.  A  comparison  of  cycles  using  He  and  N20J4  showed  a  weaker 
dependence  of  h4(<5mln)  and  Qper ( ^min)  for  than  for  He  (see 

Fig.  8). 


Fig.  8.  Dependence  of  the 
efficiency  of  a  GTP  cycle  on 
the  value  of  minimum 
temperature  head  6ml  : 

(p2  ■  10  atm;  nT  ■  0.88; 

nH  -  0.86)  I  -  T3  -  1473, 

II  -  T3  -  1273,  HI  - 

1073°K. 


The  influence  of  minimum  temperature  and  pressure  on  cycle 
effectiveness.  In  the  analysis  we  examined  the  influence  of  gas 


temperature,  having  both  auborltlcal  and  supercritical,  ahead  of  the 
compressor. 


Fig.  9.  Dependence  of  maximum 
thermal  efficiency  on  the  minimum 
and  maximum  temperatures 

and  the  pressure  after  the  com¬ 
pressor  p2:  nT  ■  0.88;  nH  ■ 

=  0.86;  0  -  0.91;  0  =  0.9; 

6min  =  10°i  T1  =  310-650°K. 

With  T1  <.  T  the  efficiency  shifts  toward  high  values  of  ir^ 
for  the  compression  ratio  in  the  compressor.  With  the  permissibly 
low  values  of  the  optimum  value  of  pressure  behind  the  compressor 
P2  is  in  the  low-pressure  region.  With  an  increase  in  T1  to  the 
range  370-430°K  the  optimum  p2  shifts  to  the  high-pressure  region 
(see  Fig.  9). 

With  T,  =  iJ50-650°K  >  T  the  optimum  again  shifts  to  the  region 

X  H  p 

of  low  values  of  p2  (see  Fig.  9).  Such  a  phenomenon  is  easily 
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explained  after  examination  of  the  behavior  of  the  effective  specific 
power  of  the  cycle  ■  f(Tj)  (Pig.  10)  and  the  work  of  gas  com¬ 
pression  in  the  compressor  AL  ■  f  (T. )  ( Fig .  11).  The  first  stage  of 

H  X 

dissociation  (N-,0^  t  2 N02)  occurs  in  the  region  of  temperatures  up  to 
450°K. 


Pig.  10.  The  dependence  of 
specific  GTP  power  on  the 

minimum  temperature  of  the  cycle 
T^  and  pressure  p2> 


Designations:  mbt/hp  -  kW/kg; 
aTa  ■  atm  (abs) . 


Fig.  11.  Dependence  of  the  work 
of  gas  compression  in  the  com¬ 
pressor  AL  on  the  minimum 

temperature  of  the  cycle  T1  and 

pressure  p2  for  a  GTP  using 

helium  and  nitrogen  tetroxide. 

Designations:  hhs/i/kt  = 

-  kcal/kg;  ata  ■  atm  (abs). 


A  rise  in  pressure  shifts  the  reaction  toward  the  region  of 
higher  temperature,  i.e.,  at  high  pressures  the  compressor  operates 
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using  a  gas  with  a  lower  gas  constant,  which  also  leads  to  a  lower 
value  of  the  compression  work  of  the  compressor. 

At  temperatures  ■  450-650°K  the  first  stage  of  the  reaction 
1 s  complete,  while  the  second  is  Just  beginning;  therefore  the 
compression  work  of  the  gas  Is  practically  Independent  of  pressure. 

The  Influence  of  the  temperature  of  the  gas  ahead  of  the  turbine 
Tj  on  the  effectiveness  of  the  cycle.  We  know  that  a  rise  in  tempera¬ 
ture  ahead  of  the  turbine  influences  the  thermal  economy  of  a  QTP, 
increasing  the  average  temperature  of  heat  input.  In  the  study,  the 
range  of  change  of  the  upper  temperature  of  the  cycle  included 
temperatures  of  600-800°C  which  are  already  used  in  power  engineering 
as  well  a3  temperatures  of  1000-1200°C  which  are  perspective  tempera¬ 
tures  . 

Analysis  showed  that  a  rise  in  gas  temperature  from  600  to  800°C 
gives  the  maximum  cycle  efficiency  increment,  approximately  10!t 
(see  Pig.  9);  with  a  further  increase  this  increment  is  significantly 
reduced : 

At]  =  i\(T,  ~  1273  *K)— i)(r,  =>  1073  °K) »  6,5V* 

An  =  n(7*,=  1473  *K)  — Tj(rg=«  l273°K)«5Vr 

This  phenomenon  can  probably  be  explained  as  follows:  with 
<  1073°K  the  second  stage  of  dissociation  2NO2  Z  2NO  +  C>2  has  its 
most  developed  nature,  i.e.,  there  Is  the  most  intense  increase  in 
the  gas  "constant"  which  is  proportional  to  the  gas  expansion  work 
in  the  turbine.  At  temperatures  >  1073°K,  however,  the  reaction 
of  dissociation  actually  terminates,  and  the  gas  "constant"  is 
stabilized  which  leads  to  a  lower  efficiency  increment.  But  at  the 
same  time,  the  work  of  regeneration  depends  on  the  pressure,  since  it 
is  in  the  region  of  the  second  stage  of  the  reaction  2N02  Z  2N0  +  02 
(Fig.  12).  Therefore,  in  this  region  low  pressure  values  are  optimum. 

As  T,  increases,  the  optimum  n.,(Ti  )  goes  toward  lower 

J.  14  K 

values  of  it  .  In  the  supercritical  region  of  temperatures  T.^  the 
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Pig.  12.  Dependence  of  gas 
expansion  work  in  the  turbine 
ALjp  on  the  compression  ratio 

it  and  pressure  p-  (T,  ■ 

-  873°K). 

Designations:  KMa/i/wr  ■ 
kcal/kg;  sts  ■  atm  (abs). 


optimum  value  it  «■  4  is  characteristic  for  the  entire  studied  pressure 

H 

range  p2  =  6-200  atm  (abs)  (Pig.  13). 


Pig.  13.  Dependence  of  psychoeffi¬ 
ciency  on  the  compression  ratio  in 
the  compressor  tr  and  the  minimum 

K 

temperature  of  the  cycle  T1  for 
pressure  p2  *  6  atm  (abs):  p2  ■ 

-  6  atm  (abs);  ■  1273°K;  nT  ■ 

■  0.88;  Hv  ■  0.86;  oo  ■  0.91; 

V„  ■  °-9°S  Atpe?  '  10°- 


The  Influence  of  the  coefficients  of  pressure  losses  on  the 
cycle  effectiveness.  Hydraulic  resistances  account  for  losses  in 
pressure  in  the  inlet  and  outlet  devices  of  the  compressor  and 
turbine  in  the  gas  line  and  the  regenerator,  in  the  cooler,  and  in 
the  heater.  Obviously  the  compression  ratio  in  the  compressor  will 
always  be  somewhat  greater  than  the  expansion  ratio  in  the  turbine; 
here  part  of  the  effective  work  of  the  cycle  is  spent  on  overcoming 
hydraulic  resistances. 


66 


The  study  was  carried  out  for  the  following  pressure  loss 
coefficient:  a  „  ■  0.85-1  along  the  low-pressure  line;  a  ■ 

H  •  B  •  A 

■  0.85-1  along  the  high-pressure  line. 

As  can  be  seen  from  Figures  14  and  15>  the  QTP  efficiency  is 
sensitive  to  pressure  lost  to  overcoming  resistances.  Consideration 
of  resistances  changes  the  optimum  compression  ratio  toward  high 
values  of  n  . 


Fig.  14. 


Fig.  14.  Dependence  of  cycle  efficiency  on 
the  coefficient  of  pressure  losses  on  the 
low-pressure  side  °H  fl:  a  —  *  l473°K;  b  — 

T3  «  1273°K;  c  -  T3  -  1073°K  (p2  -  30  atm 

(abs) ;  o  =  0.85). 

a  • 

Fig.  15.  Dependence  of  cycle  efficiency  on 
the  coefficient  of  pressure  losses  toward 
high  pressure  :  a  -  T3  =  1273°K,  ^  = 

=  330°K;  b  -  T3  =  1173°K;  c  -  T3  =  1073°K 

(T1  -  330°K,  p2  =  10  atm  (abs)). 


Analysis  of  the  results  showed  that  a  change  in  the  coefficient 
of  losses  along  the  high-pressure  line  in  the  indicated  range  has  a 
more  significant  Influence  on  the  effectiveness  of  the  cycle  (3-^) 
than  does  a  change  in  the  coefficients  of  pressure  losses  along  the 
low-pressure  line  in  this  same  range  of  change  of  o(An  =  1-1. 5/0- 
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To  compare  the  influence  of  (J  „  and  o„  „  on  the  efficiency  of 
a  QTP  using  dissociating  and  inert  (He)  working  media  we  calculated 
cycles  using  He  for  the  same  range  of  change  of  o;  these  showed  that 
the  dependence  of  cycle  efficiency  on  the  values  of  pressure  losses 
with  inert  gases  is  more  abrupt  (An  5  20%)  (Pig.  16). 


Pig.  16.  The  dependence  of 
cycle  efficiency  on  the  coeffi¬ 
cient  of  pressure  losses  o 

H  e  M 

and  o  .  in  a  QTP  using  inert 
B  •  <A 

and  chemically  reacting  working 
media:  T3  -  1073°K,  p2  - 

■  10  atm  (abs),  T1(N201j)  ■ 

■  330°K,  T1(He)  -  303°K, 


H.fl  e.fl* 

-  °-86» 


nT  -  0.88, 
-  10°. 


The  lower  influence  of  pressure  losses  along  the  line  on  the 
efficiency  of  gas  cycles  using  dissociating  working  media  allows  us 
to  permit  high  pressure  losses,  which  has  a  positive  influence  on 
the  dimensional  characteristics  of  the  device  and,  consequently, 
leads  to  better  technical-economic  characteristics  of  a  GTP  using 
dissociating  working  media  than  one  using  inert  gases. 

In  GTP's  using  dissociating  N-jO^  gas  a  decrease  in  efficiency 
of  the  turbines  has  a  much  greater  influence  on  the  effective  effi¬ 
ciency  of  the  cycle  than  in  GTP's  using  inert  gases  (Pigs.  17,  18). 
This  is  explained  by  the  high  work  ratio  of  the  gas  cycle  of  a  GTP 
using  dissociating  gas.  Higher  values  of  the  work  ratios  in  the  gas 
cycles  of  dissociating-gas  GTP's  (Pig.  19)  allows  us  to  conclude 
that  the  introduction  of  intermediate  cooling  during  compression  of 
the  gas  in  the  compressor  i3  less  effective  than  in  cycles  of  GTP's 
using  inert  gases.  Therefore,  GTP  schemes  using  intermediate  cooling 
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were  not  examined.  The  introduction,  however,  of  intermediate  super¬ 
heating  in  a  GTP  with  a  dissociating  working .medium  makes  it  possible 
to  increase  the  cycle  efficiency  by  3-5  abs.  %. 


Pig.  17. 


im  ut  tm  ut  ut  t»  t. 


Pig.  18. 


Fig.  17*  Dependence  of  cycle  efficiency  on  com¬ 
pressor  efficiency  n[(:  a  —  p2  =  10  atm  (abs); 

b  -  T3  =  1273°K;  c  -  T3  =  1073°K  (irK  =  4-7;  n  - 

=  0.9;  Tx  =  330°K). 


Fig.  18.  Dependence  of  cycle  efficiency  on  tur¬ 
bine  efficiency  nT:  a  —  T3  *  1473,  b  -  T3  =  1273, 

c  -  T3  =  1073°K;  (p2  =  10  atm  (abs);  T1  *  330°K; 

"k  -  °'78i  Smln  -  10°>- 


sw  to  mo  »  w  ns 


Fig.  19.  Dependence  of  the 
work  ratio  on  the  maximum 
temperature  for  the  cycle  of  a 
GTP  using  dissociating  and  inert 
gases . 


Our  thermodynamic  analysis  of  gas  cycles  using  a  dissociating 
working  medium  J  2N02  J  2N0  +  02  allowed  us  to  explain  their 

considerable  advantages  for  the  maximum  gas  temperatures  used  at 
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present  in  stationary  gas-turbine  construction  (600-800°C).  The  use 
of  dissociating  gases,  in  a  OTP,  in  this  temperature  range  allows  us 
to  obtain  an  effective  OTP  efficiency  of  40-45*. 

3.  A  GAS-LIQUID  CYCLE  USING 
CHEMICALLY  REACTING  SYSTEMS 

The  chemically  reacting  systems  NA**2NQ^*2NO+(\,  Al1CI^»2AIOJ. 
AljBr^naAIBr,  have  a  saturation  line  (pa-Tg)  which  is  such  that  with 
these  working  media,  simple  and  complex  gas-liquid  cycles  can  be 
constructed.  For  example,  for  the  system  NtO,**2NQ,  gas-liquid  cycle 
can  be  realized  with  the  following  minimum  parameters  in  the  cycle: 
pressure  1.1-1. 4  atm  (abs),  temperature  20-25°C  [12]  with  a  cycle 
effectiveness  which  is  3-5  abs.  %  higher  than  using  steam. 


An  outstanding  feature  of  the  gas-liquid  cycle  using  a  dissoci¬ 
ating  working  medium  compared  with  a  condensation  medium  is  the  fact 
that  it  combines  the  following  positive  qualities  of  the  thermodynamic 
cycles: 

1)  gas-turbine  -  high  mean-integral  temperature  of  heat  input; 

2)  steam-turbine  -  low  power  expended  on  pumping  the  working 
medium,  low  mean-integral  temperature  of  heat  removal; 

3)  a  cycle  using  low-boiling  working  media  -  low  heat  of  vapori¬ 
zation,  low  critical  pressure; 


4)  a  cycle  using  dissociating  working  media  -  high  effectiveness 
of  regeneration  in  the  cycle. 


These  features  allow  us  to  have  a  gas-cooled  reactor  or  a  gas- 
cooled  boiler  in  a  device  with  a  condensation  (gas-liquid)  cycle. 

The  operating  principle  of  a  device  with  a  closed  gas-liquid 
cycle  is  as  follows  (Fig.  20).  The  working  medium  in  the  liquid 
phase  enters  pump  1,  where  the  pressure  is  raised  to  the  maximum; 
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it  then  enters  regenerator-ev apor  .tor  2,  in  which  there  is  heating  of 
the  liquid,  vaporization,  and  partial  superheating  of  the  gas  due  to 
the  heat  of  the  gases  spent  in  the  turbine.  Then  the  gas  enters 
mean  heater  3»  where  it  is  superheated  to  the  maximum  temperature; 
it  then  goes  to  turbine  4.  During  the  heating  process  there  is 
dissociation  (decomposition)  of  the  initial  substance  with  the  release 
of  the  heat  of  chemical  reactions,  which  facilitates  the  isothermal 
nature  of  the  process  of  heat  feed.  The  gas,  expanding  in  the 
turbine,  accomplishes  effective  work. 


Pig.  20.  Diagram  of  a 
device  with  a  closed 
gas-liquid  cycle. 


Prom  the  turbine  the  spent  gas  enters  the  regenerator-evaporator 
and  uses  its  remaining  heat  to  heat  the  "cold"  working  medium.  The 
gas  which  is  cooled  in  the  regenerator-evaporator  enters  the  cooler- 
condenser,  where  its  condensation  and  partial  supercooling  occurs, 
and  the  cycle  is  repeated.  During  cooling  there  is  recombination  of 
the  components  into  the  initial  state  with  the  release  of  heat  of 
chemical  reaction,  which  favorably  influences  the  process  of  regen¬ 
eration  in  the  cycle. 

When  examining  possible  schemes  for  atomic  power’  plants  with 
a  fast  reactor  using  a  chemically  reacting  working  medium-coolant  we 
analyzed  the  followin  ,  schemes:  with  regeneration  using  N^O^,  with 
regeneration  and  intermediate  superheating  using  with  inter¬ 

mediate  regeneration  using  N^O^,  and  with  intermediate  regeneration 
and  intermediate  superheating  using 
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The  results  of  the  analysis  are  Given  below. 


When  analyzing  the  thfrinc  .ynamic  cycles  we  varied  the  following 
parameters:  the  pressure  p^  th«  cprr.perature  T^  at  the  turbine 
inlet;  the  pressure  at  the  turbine  nutlets;  the  turbine  efficiency 
tt—, ;  the  minimum  temperature  head  in  tne  regenerator-evaporator 

i  pal 

the  pressure  losses  along  the  low -press  ire  line  a  ;  and  the  pressure 

H 

of  intermediate  superheating. 

The  various  calculations  «■  •  .  formed  on  the  Mlnsk-22  computer 

using  the  compiled  algorithm. 

H.  ANALYSIS  OF  A  GAS-LIQUID  CYCLE 
WITH  REGENERATION 

The  Influence  of  the  pressure  and  temperature  ahead  of  the 
turbine  on  the  cycle  efficiency.  The  results  of  calculating  the 
dependence  of  cycle  efficiency  on  the  parameters  at  the  turbine  inlet 
for  the  regeneration  scheme  are  given  in  Fig.  21.  As  can  be  seen 
from  Fig.  21,  the  cycle  efficiency  increases  most  intensely  to 
pressures  of  100-1*10  atm  (abs). 


Fig.  21.  The  influence  of 
parameters  at  the  turbine  inlet 
on  the  cycle  efficiency  (with¬ 
out  intermediate  superheating): 
T:  -  290°K;  -  1.12  atm 

(abs);  -  10°;  o  -  0.8; 

°B.A  “  0,93;  nT  "  0-88;  nH  “ 

•  0.7. 

Designation:  aTa  -  atm  (abs). 


Thus,  in  the  pressure  range  10-80  atm  (abs)  with  p2  “  10  atm 
(abs)  the  cycle  efficiency  increases  by  1.5-21;  in  the  interval 
80-1*10  atm  (abs)  it  increases  by  0.5-lf;  from  a  pressure  of  1*10  atm 
(abs)  to  the  optimum  pressure  the  increase  is  ''-0.2J. 
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The  intensity  of  the  increase  in  cycle  efficiency  as  a  function 
of  temperature  at  the  turbine  inlet  and  in  the  range  =  600-1000°K 
remains  constant  and  is  ^2%  at  AT  ■  50°  (see  Pig.  21). 

In  the  temperature  range  1000-1300°K  the  increment  is  ^1?  at 
AT  =  50°. 

Such  a  change  in  cycle  efficiency  in  a  regeneration  scheme  can 
be  explained  by  the  fact  that  for  temperatures  above  1000°K  the 
chemical  reaction  ceases;  consequently,  a  change  in  regeneration 
efficiency  will  occur  only  due  to  a  change  in  heat  capacity  due  to 
temperature . 

From  an  analysis  of  the  obtained  dependences  (see  Pig.  21)  we  can 
state  that  in  schemes  having  a  cycle  with  regeneration  using  the 
chemically  reacting  mixture  t  2N02  t  2NO  +  0^  the  use  of  pressure 

greater  than  150  atm  (abs)  at  the  turbine  inlet  is  unfeasible. 


The  dependence  of  cycle  efficiency  on  counterpressure  p.. 

Figure  22  shows  the  dependence  of  efficiency,  in  a  scheme  with  regen¬ 
eration,  on  pressure  on  the  turbine  outlet.  The  intense  drop  in 
cycle  efficiency  to  a  pressure  of  5  atm  (abs)  occurs  due  to  the 
abrupt  decrease  in  operation  of  the  turbine;  with  an  increase  in 
counterpressure  from  1  to  5  atm  (abs)  the  work  of  the  turbine  decreases 
to  40  kcal/kg,  while  with  an  increase  in  counterpressure  from  5  to  10 
atm  (abs)  the  work  decreases  to  20  kcal/kg,  i.e.,  approximately 
20  kcal/wg  less  than  with  an  increase  in  counterpressure  from  1  to 
5  atm  (abs).  With  an  increase  in  counterpressure  by  more  than  10  atm 
(abs)  the  intensity  of  the  drop  in  head  in  the  turbine  is  reduced. 


Fig.  22.  Dependence  of  the  cycle 
efficiency  on  counterpressure: 


AT1111" 
per 

-  0.93;  T 


873°K;  n 
-  10° 


T 

o  ■ 

M  .A 

290°K. 


0.88;  ru  -  0.7; 


K 

0.8;  o 


B  .a 


Designation:  aTa  =  atm  (abs). 
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This  is  explained  by  the  fact  that  the  process  of  the  cycle  in 
T-S  diagram  shifts  to  the  left,  i.e.,  to  the  region  of  a  more  effec¬ 
tive  change  in  the  degree  of  dissociation  a^,  a^.  With  a  pressure 
of  1  atm  (abs)  the  degree  of  dissociation  at  the  turbine  inlet  is 
''•l,  at  the  outlet  -  'vi;  at  the  inlet  o ^  ■  0.8,  at  the  outlet  a 2  *  0.2. 
With  a  counterpressure  of  10  ^tm  (abs),  however,  the  interval  of 
change  is  smaller,  namely:  at  the  inlet  Og  =  0.8,  at  the  outlet 
a2  =  0.6  due  to  the  high  temperature  at  the  turbine  outlet.  Conse¬ 
quently,  the  average  value  of  the  gas  "constant"  in  the  turbine  in 
the  latter  case  will  be  g  ^ater,  which  favorably  influences  the  change 
of  head  in  the  turbin:. 

The  influence  -if  turbine  efficiency  on  cycle  effectiveness.  In 

the  examined  range  of  change  of  turbine  efficiency  and  parameters 
at  the  turbir?  inlet  we  obtained  a  linear  dependence  of  turbine  and 
cycle  effi  ncies.  With  a  decrease  in  turbine  efficiency  by  0.5  the 
cycle  efficiency  drops  by  'vl.Cf.  The  intensity  of  the  change  is 
independent  of  pressure. 

The  Influence  of  minimum  temperature  head  In  the  regenerator  on 

the  cycle  efficiency.  An  analysis  of  the  influence  of  minimum 

temperature  head  in  the  regenerator  (ATin*'?)  on  the  cycle  efficiency 

p»r 

is  given  in  Fig.  23.  We  obtained  a  practically  linear  dependence  of 
cycle  efficiency  on  The  Increase  in  minimum  temperature  head 

in  the  regenerator  by  10°  decreases  the  cycle  efficiency  by  0.5>. 
Selection  of  the  minimum  temperature  head  in  the  regenerator  should 
be  stipulated  by  the  weight  and  dimensional  characteristics  of  the 
regenerator-evaporator,  i.e.,  final  selection  of  ATmin 

is  uetermlned 

P  Q  l 

in  the  stage  of  technical-economic  estimates. 

On  the  basis  of  this  thermodynamic  analysis  of  gas-liquid  cycles 
with  regeneration  we  can  state  that  the  use  of  the  chemically  reacting 
mixture  N.,0^  Z  2NO2  Z  2N0  +  0^  as  the  working  medium  makes  it 
possible  to  obtain  an  efficiency  for  the  gas-liquid  cycle  with  regen¬ 
eration  at  the  pressures  and  temperatures  presently  used  in  power 
engineering  (p  <.  2h0  atm  (abs).  T  <.  838°K),  up  to  ^4  3  % . 

7H 


Pig.  23.  Influence  of  the 
minimum  temperature  head  in 
the  regenerator  on  the  cycle 
efficiency:  =  290°K;  = 

=  1.12  atm  (abs);  =  873°K; 

nT  =  0.88;  nK  =  0.7. 

Designation:  ata  =  atm  (abs). 


5.  ANALYSIS  OF  A  GAS-LIQUID  CYCLE 
WITH  REGENERATION  AND  INTERMEDIATE 
SUPERHEATING 

To  raise  the  efficiency  of  the  thermodynamic  cycle  and  the  unit 
as  a  whole  we  use  secondary  superheating  of  the  steam  (gas).  The  use 
of  intermediate  superheating  in  a  cycle  which  uses  the  chemically 
reacting  mixture  N^O^  t  2N02  t  2N0  +  02  as  the  working  medium  is  also 
expedient,  since  an  increase  of  temperature  during  expansion  increases 
the  average  gas  "constant"  in  the  turbine  due  to  the  smaller  interval 
of  change  of  the  degree  of  dissociation  and,  secondly,  increases  the 
effectiveness  of  regeneration.  The  operating  principle  of  a  cycle 
with  regeneration  and  intermediate  superheating  (Fig.  24)  is  as 
follows:  the  working  medium  in  the  liquid  phase  enters  pump  1,  Is 

"compressed"  to  the  maximum  pressure  of  the  cycle  p2  and  enters 
regenerator-evaporator  2  where  it  is  heated,  vaporized,  and  partially 
superheated  due  to  the  heat  of  the  gas  spent  in  the  turbine.  The 
superheated  working  medium  (gas)  enters  reactor  3  (intermediate  heat 
exchanger,  boiler),  where  It  is  superheated  to  the  maximum  cycle 
temperature  T^. 

The  ga3  with  maximum  temperature  T^  enters  high-pressure  turbine 
4,  is  expanded  (accomplishing  work),  and  enters  heat  exchanger 
(reactor)  3  for  secondary  superheating  of  the  gas;  this  gas  then 
enters  low-pressure  turbine  4'.  The  spent  gas  in  low-pressure  turbine 
4'  enters  regenerator-evaporator  2,  where  it  gives  its  remaining 
heat  for  heating,  evaporation,  and  partial  superheating  of  the  "cold" 
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Pig.  24.  Diagram  of  a  gas 
liquid  cycle  with  regen¬ 
eration  and  intermediate 
superheating. 


working  medium.  The  gas  which  gives  off  heat  in  regenerator-evaporator 
2  enters  cooler-condenser  5,  where  it  is  partially  recooled  to  the 
saturation  point,  condensed,  and  supercooled.  The  supercooled  working 
medium  enters  pump  1,  and  the  cycle  is  repeated. 

The  Influence  of  the  pressure  of  Intermediate  superheating  on 
cycle  efficiency.  To  perform  thermodynamic  analysis  of  a  scheme  with 
intermediate  superheating  we  must  examine  the  influence  of  the 
pressure  of  intermediate  superheating  on  the  cycle  efficiency  and 
determine  its  optimum  values,  since  the  technical-economic  calculations 
performed  using  other  working  media  (H^O,  CC^)  showed  that  the 
optimum  pressures  of  intermediate  superheating  from  thermodynamic  and 
technical-economic  calculations  coincide. 

We  determined  the  optimum  pressures  of  intermediate  superheating 
for  initial  pressures  of  80,  130,  and  240  atm  (abs)  in  the  tempera¬ 
ture  range  773-1273°K.  The  optimum  pressure  of  intermediate  super¬ 
heating,  for  all  practical  purposes,  for  all  initial  pressures  is  in 
the  range  15-30  atm  (abs).  With  a  rise  in  temperature  the  optimum 
pressure  shifts  toward  high  pressures  (Pig.  25).  At  a  temperature  of 
1173°K  and  a  pressure  of  240  atm  (abs)  the  optimum  cycle  efficiency 
with  a  pressure  intermediate  superheating  greater  than  15  atm  (abs) 
remains  practically  constant.  Such  a  flat  dependence  has  a  favorable 
influence  on  the  weight  and  dimensional  characteristics  of  the 
Installation. 
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Pig.  25.  The  dependence  of 
optimum  pressure  of  inter¬ 
mediate  superheating  p^  on  the 

temperature  of  the  gas  at  the 
turbine  inlet  with  pressure 
=  130  atm  (abs)  :  p^  = 

=  1.12  atm  (abs);  ^  =  290°K; 


nT  =  0.88;  nH  «  0.9;  °H>fl 


0.8;  0 


Vn 


0-95! 


10°. 


Designation:  ata  =  atm  (abs). 


To  estimate  the  efficiency  of  introducing  second  intermediate 
superheating  we  carried  out  various  calculations  whose  results  are 
given  in  Fig.  26.  Calculation  was  performed  for  pressure  p^  ■  2^0 
atm  (abs)  and  T^  *  838°K;  we  varied  the  pressures  of  the  first  and 
second  intermediate  superheating. 


Fig.  26.  The  influence  of  inter¬ 
mediate  superheating  pressures  on 
efficiency  with  two  intermediate 
superheatings:  p^  ■  2^0  atm 

(abs),  T,  =>  838°K. 

Designation:  aia  *  atm  (abs). 


From  an  examination  of  the  efficiency  dependence  we  see  that  the 
optimum  value  of  the  pressure  of  second  intermediate  superheat ing  p'^ 
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is  in  the  range  10-14  atm  (abs).  The  dependence  of  optimum  efficiency 
on  intermediate  superheating  pressure  p^  is  steeper  than  for  a  cycle 
with  single  intermediate  superheating;  tne  steepness  increases  with 
a  decrease  in  the  pressure  of  the  first  intermediate  superheating. 

Prom  Pig.  26  it  also  follows  that  the  values  of  the  optimum 
pressures  of  the  first  intermediate  superheating  in  a  scheme  with 
dual  intermediate  superheating  are  in  the  range  50-70  atm  (abs).  Such 
a  relationship  of  p^,  p^,  and  p|j  allows  us  to  produce  turbines  whose 
expansion  ratios  are  close  to  one  another,  for  example,  with  p^  * 

=  240  atm  (abs),  p^  *  60  atm  (abs),  p^  ■  12  atm  (abs);  ■  4,  and 

"Ten  *  5*  ■  8.5,  respectively. 

Prom  comparison  of  the  optimum  values  of  efficiency  for  single 
and  dual  intermediate  superheatings  we  see  that  in  the  latter  case 
the  cycle  efficiency  increase  by  approximately  1.4  abs.  9 . 

Prom  an  analysis  of  schemes  with  one  and  two  intermediate  super¬ 
heatings  it  follows  that  the  introduction  of  second  intermediate 
superheating  can  hardly  be  Justified  by  introducing  a  third  hot 
heat  exchanger,  which  leads  to  complication  of  the  scheme  and  might 
not  be  economically  Justified  by  a  1.4$  increase  in  thermoefficiency. 

Selection  of  second  intermediate  superheating  can  be  Justified 
only  by  technical-economic  calculations. 

The  dependence  of  cycle  efficiency  on  the  parameters  at  the 
turbine  Inlet.  Figure  27  gives  dependences  of  efficiency  on  tempera¬ 
ture  and  pressure  at  the  turbine  inlet,  from  which  it  follows  that 
the  intensity  of  a  change  of  cycle  efficiency  in  the  temperature 
range  600-1000°K  is  2  abs.  %  at  AT  -  50°.  With  a  rise  in  temperature 
the  Intensity  of  the  change  of  cycle  efficiency  drops,  and  at 
temperatures  of  10Q0-1300°K  is  1  abs.  %  at  AT  *  50°.  The  peculiarity 
of  such  a  change  is  indicated  in  an  examination  of  a  regeneration 
system  without  intermediate  superheating.  The  dependence  of  cycle 
efficiency  on  pressure  at  the  turbine  inlet  (Fig.  27)  shows  that 
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an  increase  in  pressure  ahead  of  the  turbine  above  250-260  atm  (abs) 
is  hardly  advisable,  since  the  increment  of  cycle  efficiency  at  p  ■ 
=  250-300  atm  (abs)  is  only  0.4  abs.  % . 


Pig.  27.  Influence  of  pressure 
and  temperature  at  the  turbine 
inlet  on  the  cycle  efficiency 
in  a  scheme  with  intermediate 
superheating  and  regeneration: 


p3 


240  atm  (abs);  T1  «  290°K; 

1.12  atm  (abs) ;  ATroin 

per 


«  10°;  a 

’  h.a 

=  0.93;  nT  a 

°n/n  '  °'95' 
Designation: 


0.8;  a 
0.88;  n  -  0.7; 


B  .fl 


aTa  *  atm  (abs)  . 


The  Influence  of  minimum  temperature  head  In  the  regenerator  on 
cycle  efficiency.  Figure  28  gives  the  dependence  of  cycle  efficiency 

_i  M 

on  the  minimum  temperature  head  in  the  regenerator-evaporator  ATpB|_. 

As  seen  from  Fig.  28,  we  obtain  a  practically  linear  dependence  of 
cycle  efficiency  on  AT^b”.  With  a  change  in  ATpBj?  by  10°  the  cycle 
efficiency  changes  by  1  abs.  %,  i.e.,  in  a  scheme  with  intermediate 
superheating  the  change  of  cycle  efficiency  versus  At"  is  steeper 

P  B I 

than  in  a  scheme  with  regeneration  without  intermediate  superheating 
(see  Fig.  24).  Such  a  difference  can  be  explained  by  the  greater 
efficiency  of  regeneration  in  a  cycle  with  intermediate  superheating, 
since  the  introduction  of  Intermediate  superheating  raises  the 
temperature  at  the  regenerator  inlet  on  the  low-pressure  side,  which 
leads  to  an  increase  in  regeneration  due  to  the  difference  QAa  for 
high  and  low  pressures.  With  a  change  in  minimum  head,  however,  the 
effect  of  a  change  In  the  diflerence  of  the  heats  of  chemical  reactions 
will  vary.  In  a  scheme  without  intermediate  superheating  it  will  be 
less,  since  the  temperature  at  the  regenerator  inlet  is  lower,  while 
che  temperature  on  the  low-pressure  side  at  the  regenerator  outlet 
is  Identical  for  both  schemes. 


I 
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Pig.  28.  The  influence  of  minimum 
temperature  head  in  the  regenerator  on 
the  cycle  efficiency:  1  —  ■  973°K, 

P3  -  130  atm  (abs);  2  -  T3  =  838°K,  p3  = 

-  2l»0  atm  (abs);  3  -  T3  =  873°K,  p3  = 

■  430  atm  (abs)  (with  p3  =  20  atm  (abs), 
nT  -  0.88;  nK  -  0.7;  CH.a  ■  0.8;  Pl  . 

=  1.12  atm  (abs);  T.  =  290°K;  a 

I  B  .fl 

■  0.93). 


The  influence  of  counterpressure  p^  on  cycle  efficiency.  The 
efficiency  of  a  cycle  depends  on  the  compression  ratio  in  the  cycle 
as  well  as  the  absolute  values  of  the  counterpressure.  Figure  29 
shows  the  dependence  of  cycle  efficiency  on  counterpressure  p^.  As 
can  be  seen  from  Fig.  29  the  cycle  efficiency  depends  on  the  change 
in  counterpressure.  Therefore,  selection  of  counterpressure  should 
be  determined  by  the  purpose  of  the  installation,  since,  first,  the 
Increase  in  pressure  behind  the  turbine  leads  to  a  reduction  in  size 
of  the  regenerator-evaporator,  and  also  to  a  decrease  in  the  exhaust 
area  and  the  diameter  of  the  exhaust  lines,  and  second,  an  increase 
in  counterpressure  makes  it  possible  to  use  an  economic  quantitative 
method  for  controlling  the  installation  in  partial-power  regimes. 


Fig.  29.  Dependence  of  cycle  effi¬ 
ciency  on  counterpressure:  1  —  T3  ■ 


973°K,  p, 


130  atm  (abs) ;  2  —  T. 


838°K,  p3  ■  240  atm  (abs);  3  -  T3 
873°K,  p.j  ■  130  atm  (abs)  (with  nT 

■  0. 8;  a  - 

20  atm  (abs) 


-  0.88;  n  -  0.7;  0 


H 


95;  p^ 


■  °-93-  4n/n 
Wpe?  ■  1°°>- 
Designation:  ata  ■  atm  (abs). 


The  minimum  value  of  counterpressure  is  determined  by  the 
temperature  of  the  cooling  medium. 
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Figure  30  gives  the  dependence  of  cycle  efficiency  on  turbine 
efficiency.  In  giving  the  dependence  of  cycle  efficiency  on  nT 
assumed  that  hjg^  =  nTH^.  We  obtained  a  practically  linear  dependence. 
The  nature  of  change  is  practically  independent  of  temperature  and 
pressure  at  the  turbine  inlet. 


Fig.  30.  The  influence  of  turbine 
efficiency  on  cycle  efficiency:  1  — 

T3  =  973°K,  p3  “  130  atm  (abs);  2  - 

T3  ■  838°K,  p3  *  240  atm  (abs);  3  — 

T3  -  873°K,  p3  *  130  atm  (abs)  (px  - 

1.12  atm  (abs);  T,  -  290°K;  n  -  0.8; 

i.  H  • 


V*  ■  “-93;  -  10° i 

*  0.1;  p^  ■  20  atm  (abs); 


nA;  ne  - 
-  0.95). 


Figure  31  gives  the  results  of  calculating  the  influence  of 
pressure  losses  in  the  low-pressure  line  on  cycle  efficiency.  From 
Fig.  31  it  follows  that  with  use  of  the  chemically  reacting  mixture 
t  2N02  t  2N0  +  02  as  the  working  medium  of  the  cycle,  pressure 
losses  up  to  20?  reduce  the  cycle  efficiency  by  0.5  abs.  ?. 


Fig.  31.  Dependence  of  cycle 
efficiency  on  pressure  losses: 

1  -  T3  -  973°K,  p3  -  130  atm  (abs); 

2  —  T3  •  838°K,  p3  -  240  atm  (abs); 

3  -  T3  -  873°K,  p3  -  130  atm  (abs) 

(p^  ■  20  atm  (abs);  p^  »  1.12; 

T 1  -  290°K,  nT  -  0.88;  nK  -  0.7; 


rftln  m 
per 

-  0.95) 


10c 


B  .fl 


0,93;  °n/n 


On  the  basis  of  thermodynamic  analysis  of  cycle  with  single  and 
dual  intermediate  heating  we  can  draw  the  following  conclusions: 


1)  with  an  increase  in  pressure  the  relative  increase  in  effi¬ 
ciency  and  specific  work  Increases,  which  limits  the  feasibility  of 
using  high  pressures; 
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2)  with  a  rise  in  temperature  ahead  of  the  turbine,  the  cycle 
efficiency  sharply  increases.  An  increase  in  temperature  by  100° 

(in  the  range  =  600-1200°K)  raises  the  cycle  efficiency  by  4  abs. 
%,  regardless  of  the  pressure.  At  temperatures  above  1200°K  the 
intensity  of  an  increase  in  efficiency  decreases,  and  is  1.5-2?. 
Consequently ,  it  is  expedient  to  use  the  chemically  reacting  mixture 

n2°4  *  2N02  *  2N0  +  °2  as  a  working  medium  for  temperatures  of  600- 
1200°K . 

The  disadvantages  of  these  cycles  using  N^  include  the  rela¬ 
tively  high  temperature  of  the  gas  at  the  inlet  to  the  cooler- 
condenser  (about  80°C).  In  the  recooler  we  get  a  heat  of  approxi¬ 
mately  90 %  of  the  condensation  heat.  To  reduce  this  temperature,  and 
also  to  reduce  the  dimensions  of  the  regenerator-evaporator  we 
analyzed  cycles  with  intermediate  regeneration. 

6.  A  GAS-LIQUID  CYCLE  WITH  INTERMEDIATE 
REGENERATION 

Plgure  32  shows  a  diagram  of  a  cycle  with  intermediate  regener¬ 
ation.  A  feature  of  this  cycle  is  that  in  the  high-pressure  turbine 
(TBA)  [H-PT]  the  gas  is  expanded  to  a  certain  intermediate  pressure, 
after  which  it  enters  the  regenerator.  In  the  regenerator  the  gas 
gives  off  part  of  its  heat  for  heating,  evaporation,  and  superheating 
of  the  coolant  coming  from  the  pump;  after  this  it  goes  to  the  low- 
pressure  turbine  (THU)  [L-PT]  where  it  is  expanded  to  the  minimum 
pressure  of  the  cycle  and  enters  the  condenser.  A  rise  in  pressure 
and  temperature  of  the  gas  on  the  hot  side  of  the  regenerator  is 
reflected  In  an  Improvement  of  the  weight  and  dimensional  character¬ 
istics  of  the  regenerator  and  the  regeneration  process. 

i  he  pressure  of  the  Intermediate  regeneration  was  selected 
based  on  the  following  conditions. 

First,  with  an  increase  in  pressure  at  the  H-PT  outlet  there  is 
an  Increase  in  enthalpy  of  the  gas  at  outlet  from  the  H-PT,  which 
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Pig.  32.  T-S  diagram  of  a 
cycle  with  intermediate 
regeneration. 


increases  the  head  in  the  regenerator;  this  increase  in  head  is  more 
intense  than  the  increase  in  enthalpy  of  the  gas  at  the  outlet  from 
the  H-PT.  Thus,  the  enthalpy  of  at  the  outlet  from  the  regen¬ 

erator  on  the  hot  side  decreases.  Prom  the  t-H  diagram  it  follows 
that  a  decrease  in  enthalpy  of  the  gas  at  the  regenerator  outlet  on 
the  hot  side  together  with  a  rise  in  pressure  results  in  approach  to 
the  saturation  line,  i.e.,  at  the  inlet  to  the  L-PT  turbine  saturated 
vapor  may  appear. 

Second,  increasing  the  pressure  in  the  regenerator  on  the  hot 
side,  we  thus  lower  the  enthalpy  of  the  gas  at  the  point  of  minimum 
temperature  head  in  the  regenerator  on  the  hot  side  H^p,  which  leads 
to  an  increase  in  regeneration  heat. 

Third,  with  an  increase  on  the  hot  side  of  the  regenerator  the 
steepness  of  the  isobar  in  the  t-H  diagram  decreases,  i.e.,  the 
point  of  minimum  temperature  head  in  the  regenerator  shifts  to  the 
region  of  high  temperatures;  consequently,  the  enthalpy  of  at 

point  ATmln  0n  the  cold  side  Hx  increases,  which  increases  the  heat 
of  regeneration. 

Expansion  in  the  L-PT  occurs  in  the  region  of  isobars  close  to 
one  another;  therefore  the  total  work  of  the  turbine  decreases  com¬ 
pared  with  that  during  ordinary  regeneration.  We  examined  schemes 
with  regeneration  and  intermediate  superheating  and  with  regeneration 
without  intermediate  superheating. 
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The  maximum  increase  in  cycle  efficiency  without  intermediate 
superheating  and  with  intermediate  regeneration  was  3.5  abs.  t  and 
only  1.2  abs.  %  lower  than  in  a  cycle  with  intermediate  superheating 
without  intermediate  regeneration.  The  optimum  pressure  of  inter¬ 
mediate  regeneration  in  the  given  scheme  with  initial  pressure  240  atm 
(abs)  pnp  =  15-20  atm  (abs);  with  an  initial  pressure  of  130  atm  (abs) 
it  was  10-15  atm  (abs). 

The  introduction  of  intermediate  regeneration  into  the  scheme  of 
a  cycle  with  intermediate  superheating  also  causes  a  rise  in  cycle 
efficiency,  approximately  1  abs.  X,  but  in  this  case  the  optimum 
pressure  of  intermediate  regeneration  reduces  to  8-12  atm  (abs). 

To  analyze  the  efficiency  of  a  cycle  with  intermediate  regener¬ 
ation  we  also  examined  the  influence  of  pressure  losses  in  the 
regenerator  on  cycle  efficiency.  Figure  33  gives  the  results  of 
these  calculations  and  compares  them  with  the  efficiency  of  an 
ordinary  cycle. 


Fig.  33»  The  influence  of 
pressure  losses  in  the  regener¬ 
ator  on  the  hot  side  on  cycle 
efficiency . 

Designation:  aTa  ■  atm  (abs). 


As  can  be  seen  from  this  dependence  (Fig.  3*0  >  the  efficiency 
of  cycles  with  intermediate  regeneration  is  higher  than  those  using 
steam,  in  the  entire  investigated  temperature  range. 

These  investigations  of  thermodynamic  cycles  usin^  chemically 
reacting  working  media  show  that  their  use  opens  new  possibilities 
for  improving  the  technical-economic  indices  of  power  plants. 


♦ 
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Pig.  Comparison  of  the 

ef ficiency • of  cycles  using  the 
chemically  reacting  mixture 
N201|  t  2N02  t  2N0  +  02  and 

steam:  1  —  cycle  with  inter¬ 

mediate  superheating  and  inter¬ 
mediate  regeneration  using 
N20^,  2  —  cycle  with  inter¬ 
mediate  superheating  using  N2Oi< 

3  -  cycle  with  intermediate 
superheating  using  H20,  b  — 

simple  cycle  with  regeneration 
using  N20^  (p^  =  2^0  atm  (abs)) 
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APPENDIX 


Table  1.  Density,  enthalpy,  and  entropy 
in  the  saturation  line  ?  2NO2 


t.  "C 

9.  om 

p*\ 

fr,i 

SM. 

kxcu/xx- 

"r. 

sr. 

fJUUlU- 

iJCM • 

t/CM 1 

tUUU/Ki 

•tpad 

JUUM/M 

• V •» 

0,00473 

0,00571 

0,00692 

0,00o32 

0,01000 

0,01198 

0,01422 

0,01688 

0,01996 

0,02329 

0,02760 

0,03231 

0,03780 

0,04430 

0,05140 

0,05955 

0,06890 

0,08025 

0,09325 

0,10930 

0,12/00 

0,14820 

0,17520 

0,21140 

0,26200 

0.35500 

0,55000 


3 

Designations:  an  =  atm  (tech);  r/cm  ■ 

=  g/cnr;  Kwa/i/nr  =  kcal/kg;  KKa/i/Kr*rpa,q 
=  kcal/kg*deg. 


Table  2.  Enthalpy,  entropy,  heat  capacity, 
and  density  in  the  gas  phase.  Designations: 
KKa/i/ Hr  =  kcal/kg;  KKa/i/urTpafl  »  kcal/kg  x 
x  deg;  wr/cn  =  kg/cm. 


u  *c 

H,  — 

KM 

s. 

KM'lpQd 

30 

157,37 

P  —  1  K SlCjfl 

0,9142 

40 

172,91 

0,9647 

50 

191,06 

1,0217 

60 

210,62 

1,0813 

70 

229,72 

1,1378 

80 

246,53 

1,1862 

90 

260,08 

1,2240 

100 

270,39 

1,2521 

110 

278,06 

1,2724 

120v 

285,08 

1,2907 

130 

*  290,00 

1,3031 

140 

294,23 

1,3134 
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Table  2  (continued). 


1.  *c 

H 

s  .  "" 

uu 

CM* 

u 

utpcd 

r  u<4W 

900 

774,94 

1,9096 

0,33808 

1,248 

950 

791 .43 

1,9231 

0,32264 

1,194 

1000 

807,28 

1,9361 

0,31205 

1,144 

1050 

822,69 

1,9180 

0,30474 

1,099 

1100 

837,80 

1,9592 

0,29967 

1,057 

1150 

852,69 

1,9698 

0.29GI6 

1,019 

1200 

867,43 

1,9800 

0,29373 

0,983 

p  =  S  m/ck* 

9 

70 

182,97 

0,9502 

1,58890 

12,290 

80 

199,55 

0,9979 

1,70180 

10,990 

90 

216,79 

1,0161 

1,70970 

9,897 

100 

233,49 

1,0916 

1,59630 

9,005 

no 

248,59 

1,1316 

1,39160 

8,294 

120 

262,14 

1,1666 

1,17370 

7,726 

130 

272,82 

1,1935 

0,91494 

7,294 

MO 

281,37 

1,2145 

0,75826 

6,950 

150 

288,27 

! ,231 1 

0,61876 

6,670 

160 

293,97 

1,2444 

0,52287 

6,436 

170 

298,87 

1,2556 

0,46021 

6,232 

180 

303,27 

1,2654 

0,42184 

6,054 

190 

307,39 

1,2744 

0,40196 

5,890 

200 

311,37 

1,2829 

0,39179 

5,739 

250 

332,34 

1,3250 

0,46114 

5,087 

300 

359,16 

1,3739 

0,61583 

4,523 

350 

394,41 

1,4329 

0,79133 

4,008 

400 

438,14 

1,5003 

0,94588 

3,538 

450 

487,68 

1,5712 

1,01920 

3,124 

500 

533,35 

1,6390 

0,99266 

2,775 

550 

585,60 

1,6982 

0,88888 

2,493 

600 

626,73 

1,7467 

0,75592 

2,269 

650 

661,30 

1,7852 

.0,63107 

2,092 

700 

690,22 

1,8158 

0,53101 

1,948 

750 

714,82 

1,8104 

0,45733 

1,829 

800 

736,30 

1,8610 

0,40539 

1,727 

850 

755,61 

1 ,8786 

0,36916 

1,639 

900 

773,43 

1,8911 

0,34476 

1,562 

950 

790,21 

1,9081 

0,32775 

1,493 

1000 

806,29 

1,9210 

0,31598 

1,430 

1050 

821 .87 

1,9330 

0.30760 

1,374 

1100 

837,11 

1,9143 

0,30203 

1,322 

1150 

852,11 

1,9550 

0,29807 

1,274 

1200 

866,94 

1,9653 

0.29527 

1,229 

p  =  6  Kite*1 

70 

178,06 

0,9309 

1,51240 

15,320 

80 

194,00 

0,9768 

1,65420 

13,660 

90 

220,91 

1,0240 

1,69210 

12,260 

100 

227,71 

1,0698 

1,62790 

11,110 

110 

243,33 

1,1112 

1,46080 

10,190 

120 

257.65 

1,1482 

1,25810 

9.452 

130 

269,20 

1,1773 

1,02910 

8,885 

140 

278.58 

1,2003 

0,82954 

8,435 

150 

286,13 

1,2184 

0,67432 

8,072 

160 

292,32 

1,2329 

0,56382 

7,772 

170 

297,58 

1.2450 

0,48905 

7,515 

180 

302,23 

1,2553 

0,44100 

7,292 

190 

306,50 

1,2647 

0,41384 

7,090 

200 

310,57 

1,2734 

0,40104 

6,905 

250 

331,41 

1,3152 

0,45468 

6,115 

300 

357,51 

1,3628 

0,59685 

5,441 

350 

391,64 

1,4199 

.  0,76887 

4,829 

400 

434,07 

1,4853 

0,92000 

4,269 

450 

482,49 

1,5546. 

1,00)4 

3,774 

500 

532,61 

1,6216 

0.9S847 

3,353 

550 

579,98 

1,6810 

0,89766 

3,010 

600 

621.72 

1,7303 

0,77111 

2,738 

650 

657,10 

1,7697 

0,61786 

2,521 

90 
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Table  2  (continued). 


/.  *c 

H,  ^ 

KS 

S.  — 

Kt-tpai 

'  urp t 

P.  — 

CM 

GOO 

613.52 

1,7041 

0.79255 

• 

3,682 

G50 

650,12 

1,7448 

0.67348 

3,384 

700 

CHI.lt 

1.7776 

0,57058 

3,143 

750 

707,54 

1,8041 

0,49000 

2,944 

600 

730.51 

1 .8200 

0,43104 

2,776 

850 

751,00 

1 .8447 

0,39014 

2,631 

900 

769,73 

1,8610 

0,36072 

2,504 

950 

787,22 

1,8756 

0,34007 

2,391 

1000 

803,84 

1,8869 

0,32554 

2,290 

1050 

819,85 

1,9013 

0,31527 

2,198 

1100 

835,42 

1,9128 

0,30798 

2,114 

1150 

850,68 

1,9237 

0,30277 

2.037 

1200 

865,72 

1,9341 

0,29905 

1.965 

p  =  9  m/m* 


181,86 

0,9310 

1,50480 

197,67 

0,9753 

1,61670 

214,17 

1,0201 

1,63830 

230,39 

1,0631 

1,55920 

246,05 

1,1036 

1,42200 

250,50 

1,1375 

1,21320 

277.82 

1,1653 

1,00470 

280,07 

1 , 1875 

0,82369 

287.61 

1,2052 

0,67891 

293,90 

1,2196 

0,57409 

299,28 

1,2316 

0,50071 

.  304,05 

1,2421 

.  0,45381 

308,46 

1,2515 

0,42560 

329,37 

1,2935 

0.43831 

354,09 

1,3385 

0,55859 

385.91 

1,3918 

0,71555 

425.53 

1,4529 

0,86295 

471,39 

1 ,5185 

0,95789 

.519.98  - 

1,5835 

0,97197 

567,26 

1.6427 

0,90902 

610,07 

1,6933 

0,80032 

647,13 

'  1,7346 

0,68355 

678,63 

1,7678 

0,58055 

705,53 

1,7948 

0,49930 

728,90 

1 .8171 

0,43906 

749,70 

1,8361 

0,39578 

768,68 

1,8526 

0,36513 

786,37 

1 ,8674 

0,34350 

803.14 

1,8808 

0,32822 

819.27 

1,8933 

0,31738 

834,93 

1,9049 

0.30965 

850,27 

1,9159 

0,30411 

865,38 

1.9263 

p  —  10  «/«*• 

0,30014 

80 

178.73 

0,9194 

1,47470 

90 

194.22 

0,9627 

1,58480 

100 

210,52 

1,0071 

1,62540 

no 

226.78 

1,0502 

1,57000 

120 

242.66 

1,0912 

1,45380 

130 

256, GO 

1,1263 

1,25860 

140 

268,41 

1,1554 

1,05210 

150 

278,14 

1,1788 

0,86719 

160 

286,10 

1,1974 

0,71446 

170 

292,71 

1,2125 

0,60162 

180 

298,33 

1,2251 

0,52070 

190 

303,29 

1,2359 

0,46787 

200 

307,  80 

1,2456 

0,43487 

250 

336,83 

1,2878 

0,43520 

300 

353,25 

1,3323 

0,54955 

L^T^. 
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Table  2  (continued). 


• 

1.  «c 

wu 

a 

KJ 

r 

350 

384,51 

1,3846 

400 

423,42 

1,4446 

450 

468,61 

1,5093 

500 

516,76 

1,5737 

550 

563,92 

1,6328 

600 

606,95 

1,6836 

650 

644,39 

1,7253 

700 

676,35 

1,7590 

750 

703,66 

1,7864 

800 

<27,40 

1,8091 

850 

748,49 

1,8283 

900 

767,70 

1,8451 

950 

785,57 

1,8600 

1000 

802,48 

1,8736 

1050 

818,72 

1,8861 

1100 

834,48 

1 ,8978 

1150 

849.89 

1,9088 

1200 

865,05 

1,9193 

0,70247 

0,81836 

0,94093 

0,90618 

0.91035 

0,80674 

0,69230 

0,58943 

0,50718 

0,44558 

0,40099 

0,36923 

0,34671 

0,33073 

0,31936 

0,31122 

0,30537 

0,30115 


110 

120 

130 

140 

150 

160 

170 

180 

190 

200 

250 

300 

350 

400 

450 

500 

550 

600 

650 

700 

750 

800 

850 

900 

950 

1000 

1050 

1100 

1150 

1200 


120 

130 

140 

150 

160 

170 

180 

190 

200 

250 

300 

350 

400 

450 

500 

550 


200.71 
216,96 
232,67 
247,28 

200.32 
271,60 
281,14 
289,16 
295.94 
301.78 

324.88 

347.89 
376,02 

410.72 
451.51 

496.33 
542,13 
585,87 
625,50 
o60,27 
690.43 
716,69 
739,91 
700,82 
7S0.04 
798,01 
815,08 
831,48 
847.39 
862,93 


p  =  20  a/CJfi 

i  0,9607 
1,0027 
1,0123 
1 ,0783 
1,1096 
1,1362 
1,1581 
1,1761 
1,1910 
1,2035 
1,2500 
1,2919 
1,3389 
1,3925 
1,4509 
1,5103 
1,5682 
1,6199 
1.6641 
1,7008 
1,7311 
1,7562 
1,7775 
1 ,7957 
1,8118 
1,8263 
1,8395 
1,8516 
1,8630 
1,8737 
p= 30 


200,61 

216,49 

231,80 

246,16 

259,12 

270.52 
280,35 

288.72 
295,87 
321,89 

344.72 
.371,42 

403,99 

442,39 

485,08 

529.53 


0,9494 

0,9893 

1,0270 

1,0616 

1,0921 

1,1183 

1,1404 

1,1583 

1,1742 

1,2268 

1,2685 

1,3130 

1,3633 

1,4183 

1,4754 

1,5311 


1,56360 

1,56150 

1,47990 

1,34700 

1,18140 

1,00630 

0,84585. 

0,71281 

0,61080 

0,53737 

0,43053 

0,50127 

0,62567 

0,75815 

0,86428 

0,91562 

0,90302 

0,83684 

0,74304 

0,64595 

0,56021 

0,49117 

0,43835 

0.39910 

0,37032 

0,34930 

0,33394 

0,32265 

0,31430 

0.30807 


1,56140 
1,54160 
1,45100 
1,32950 
1,18230 
1,02430 
1 ,87445 
0,74628 
0,64418 
0,44525 
0,48205 
0,58830 
0,71044 
0,81647 
0,87869 
0,86681 


8,126 

7,213 

6,400 

5,695 

5,108 

4,634 

4,253 

3,945 

3,691 

3,477 

3,293 

3,133 

2,990 

2,863 

2,747 

2,641 

2,544 

2,455 


44,130 

39,350 

35,590 

32.650 
30,280 

28.410 
26,890 

25.650 

24.630 
23,780 
20,680 
18,390 

16.410 

14.630 
13,040 
11,640 
10,440 

9,453 

8.645 

7.987 
7,447 
6,996 
6.612 
6,280 

5.988 
5.727 
5,492 
5,278 
5.082 
4,901 


67,880 

60,010 

53,970 

49,270 

45,510 

42.560 
40,110 
38,150 

36.560 
31,280 
27,710 

24.730 
22,080 

13.730 
17,640 
15,840 
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Table  2  (continued). 


H  **" 

S  -  **". 

Wi 

« 

nj-ipmd 

'■  „  «prf 

lor.o 

806,92 

1,7738 

0,30095 

13,690 

1100 

824,56 

1,7870 

0,34467 

13,150 

1150 

841,49 

1,7991 

0,33237 

12,650 

1200 

857,87 

1,8105 

0,32305 

12,180 

150 

205,47 

p- 60  ta/cu* 
0.9438 

1,50290 

136,000 

ICO 

221,18 

0.9809 

1,43320 

117,300 

170 

230,52 

1,0165 

1,37.360 

104,100 

180 

250.94 

1,0194 

1,24850 

94, .370 

ISO 

263,96 

1,0784 

1,10630 

86,800 

200 

275,46 

1,1034 

0.96718 

80,770 

250 

313,91 

1,1821 

0,52782 

64,130 

300 

338,59 

1,2272 

0,47282 

55,870 

350 

363,83 

1,2693 

0,53893 

49,680 

400 

393,35 

1.3148 

0.63865 

44,360 

450 

427,94 

1,3644 

0,73785 

39,740 

500 

466,88 

1,4165 

0,60857 

35,660 

550 

508,43 

1,4686 

0,84114 

32,040 

600 

550,57 

1,5182 

0,83146 

28,970 

650 

591 ,21 

1,5635 

0,78238 

26,430 

700 

628.83 

1,6033 

0,71350 

24,340 

750 

662,75 

1 ,R373 

0,63772 

22,640 

800 

692,95 

1,6663 

0,56679 

21.230 

850 

719,84 

1,6909 

0,50605 

20,040 

900 

743,95 

1,7120 

0,45684 

19,030 

950 

765,86 

1,7304 

0,41827 

18,140 

1000 

786,05 

1,7467 

0,38865 

17,360 

1050v 

804,93 

1 .7613 

0,36615 

16,640 

1100 

*  822,81 

1,7747 

0,34915 

15,990 

1150 

839,94 

1,7870 

0,33634 

15,370 

1200  • 

856,50 

1.7985 

0,32671 

14,800 

150 

191,31 

p- 70  Kija i* 
0,9036 

1,68460 

190.000 

160 

•  208,06  - 

0,9460 

1,46780 

155,100 

170- 

224,07 

0,9833 

1,41530 

133,500 

180 

239,62 

1,0187 

1,33900 

118,300 

190 

254,17 

-  1,0514 

1,21180 

107,100 

200 

267,24 

1,0799 

1,07120 

98,420 

250 

311,14 

1,1702 

0,56513 

75,650 

300 

377,01 

1,2176 

0,47930 

65,380 

350 

362,24 

1,2597 

0,53205 

57,980 

400 

391,27 

1,3045 

0,62518 

51 ,750 

450 

425,16 

1,3531 

0,72226 

46,370 

500 

463,40 

1,4042 

0,79493 

41,620 

550 

504,41 

1,4556 

0,83132 

37,380 

COO 

546,25 

1,5050 

0,82629 

33,800 

650 

586.86 

1,5503 

0,78297 

30,820 

700 

624,66 

1,5903 

0,71668 

28,380 

750 

658.88 

1 ,6247 

0,64270 

26,390 

800 

689.43 

1,6541 

0,57224 

24,730 

850 

716,64 

1,6791 

0,51118 

93,330 

900 

741.04 

1,7006 

0,46134 

>2,160 

950 

763,18 

1,7192 

0,42215 

21,130 

1000 

783,57 

1,7358 

0,39200 

20,210 

1050 

802,61 

1,7506 

0,36930 

.19,370 

1100 

820,64 

1,7641 

0,35226 

18,610 

1150 

837,92 

1.7765 

0,33961 

17,890 

1200 

854.64 

1,7880 

0,33038 

17,210 

• 

150 

176,85 

P---80 

0,8699 

2,43760 

288,4(10 

160 

196,93 

0,9170 

1,60180 

203,600 

170 

213,79 

0,9.560 

1,46400 

167,800 

180 

229,96 

0,9929 

1,38310 

145,200 

95 
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Table  2  (continued). 


1.  *c 

If  «*" 

s  ““ 

(Ski 

,  ■ 

ftf 

’  mt-tpti 

€M 

-  800 

680,73 

1 ,6242 

0,59009 

35,250 

850 

708,89  - 

1,6499 

0,52913 

33,310 

900 

734,19 

1,6720 

0,47792 

31,660 

950 

757,13 

1,6912 

0,43680 

30,220 

1000 

778,19 

.  1.7081 

0,40474 

28,940 

1050 

797.81 

1,7233 

0,38027 

27,770 

1100 

816,33 

1,7371 

0,36197 

26,660 

1150 

834,03 

1 ,7499 

0,34838 

25,650 

1200 

851,12 

1,7617 

0,33852 

24,660 

FTD-MT-24- 192^-71 


p=  150  m/cM* 


143,78 

0,78675 

1,58050 

161,99 

0,82750 

1,79810 

182/0 

0,87246 

1,91460 

201.  ,i 

0,91412 

1,61060 

218,63 

0,95125 

1,45380 

286,62 

1,0920 

0,89(44 

325,18 

1,1637 

0,57658 

.352,97 

1,2102 

0,52419 

380,58 

1,2527 

0,57096 

411, 62 

1,2971 

0,65395 

446.55 

1,3440 

0,72474 

484,55 

1,3916 

0,77370 

524,31 

1,4384 

0,79275 

561,28 

1,4829 

0,77633 

602.75 

1,5236 

0,73223 

633,54 

1,5596 

0,67221 

671,10 

1,5909 

0,60776 

700,40 

1,6179 

0,54707 

726,74 

1 ,6412 

0,49449 

750,58 

1,6615 

0,45137 

772,40 

1,6793 

0,41735 

792,62 

1,6951 

0,39125 

811,64 

1,7095 

0,37214 

829,76 

1,7226 

0,35778 

847,21 

1,7347 

0,34758 

p=200  title* 

132,74 

0,7591 

1,25950 

(46,65 

0,7898 

*  - 100 

162,50 

0,8244 

1  J 

179,93 

0,8620 

1  .uO 

197,81 

0,9001 

1 ,55610 

272,47 

1,0536 

1,05010 

318,18 

1,1391 

0,63919 

348,11 

1,1895 

0,53902 

376, 15 

1,2327 

0,56809 

406,62 

1,2762 

0,63232 

440,52 

1,3215 

0,69904 

477,67 

1,3680 

0,75336 

516,88 

1,4141 

0,77734 

556,66 

1,4584 

0,76917 

595,37 

1,4993 

0,73337 

631,76 

1,5359 

0,67972 

665,14 

1,5081 

0,61896 

695,32 

1,5960 

0,55963 

722,49 

1,6202 

0,50687 

747,04 

1,6411 

0,46281 

769,43 

1,6595 

0,42758  . 

790,11 

1 ,6759 

0,40033 

809,46 

1,6905 

0,37980 

827,79 

1,7038 

0,. 36475 

815,39 

1,7160 

0,35411 

p— 250  title* 

127,32 

0,7443 

1,06600 

139,33 

0,7708 

1,19770 

152,60 

0,7998 

1,31500 

167,62 

0,8317 

1,46780 

183,50 

0,8653 

1,48750 

259,87 

1,0211. 

1,14400 

97 

786.600 
631,000 
459,000 

356.500 

296.600 
178,000 

142,300 

123.100 

109.100 
97.630 
87,760 
79.020 
71,510 
65.190 
59.950 
55,640 
52.090 
19.130 
46.610 
44.430 

42.500 
40,750 
39,120 

37.600 
36,140 

892.900 

802.900 

692.200 

572.100 

469.100 

251.900 

192.200 

163.100 

144.100 

128,700 
115.400 
104,000 
94,230 
85,840 
78,850 
73,070 
68,301 
64,330 
GO, 980 
58,080 
55,530 
53,220 
51,090 

49.100 

47.200 


944.300 

884.700 

802.700 

716.300 
624,800 

331.200 


Table  2  (continued). 


».  *c 

u 

r  KJUU 

KttftJ 

P.  — 

CM 

« 

Ki-ipad 

**  Ki  ipad 

300 

310.65 

1,1164 

0,73124 

211,500 

CO 

3M.M 

1,1727 

0,36168 

202,600 

400 

372.06 

1,2167 

0,30016 

177,800 

ro 

402,58 

1,2595 

0,01864 

158,(730 

0X3 

4.15,05 

1,3637 

0.08I0G 

142,400 

550 

471 .64 

1,3190 

0,7. 1573 

128,400 

600 

510.19 

1,3940 

0,76312 

116,300 

630 

549,44 

1,4376 

0,76150 

106,000 

700 

56b  ,  00 

1,4764 

0,7,1232 

97,380 

730 

024,36 

1,5153 

0,68.399 

90,280 

800 

60S, 33 

1,5480 

0,62600 

81,450 

830 

669,01 

1,5766 

0,56877 

79,630 

OOO 

716,72 

1,6015 

0,51619 

75,560 

930 

741,79 

1 ,6232 

0,47158 

72,070 

1000 

761,64 

1,6122 

0,43556 

68,980 

1050 

765,72 

1,6590 

0,40758 

66.170 

MOO 

605,43 

1,6741 

0,38055 

63,570 

1150 

824,09 

1,6878 

0,37130 

61,090 

1200 

841,99 

1,7002 

0,36068 

58,700 

P-300  Ki/CM* 

IGO 

123,92 

0,7340 

1 ,96242 

931,300 

170 

134,06 

0,7578 

1,05090 

929,300 

IH3V 

146  64 

0.7840 

1,17810 

869,600 

100 

•  ijy,<Jl 

0,8120 

1,32060 

795,600 

200 

174,48 

0,6126 

1,35070 

722,800 

250 

249,13 

0,9941 

1,19510 

410,300 

300 

303,96 

1,0965 

0,80864 

289,500 

350 

310, '26 

1,1579 

0,59145 

240,600 

400 

369.41 

1,2029 

0,55628 

210,400 

450 

399,02 

1,2153 

0,61033 

187,500 

500 

431,48  “ 

1,2687 

0,66757 

168,500 

530- 

466,87 

1,3329 

0.72214 

152,100 

COO 

501,35 

1,3770 

0,74944 

138,000 

650 

512,71 

'  1.4195 

0,74973 

125,900 

700 

560,49 

1 ,4593 

0,72425 

115,700 

730 

616,48 

1,4954 

0,68004 

107,300 

600 

649,90 

1,5274 

0,62622 

100,300 

650 

660,44 

1,5555 

0,57098 

94,580 

900 

703,10 

1,5800 

0.52004 

69,690 

930 

7.33,34 

1,0015 

0,17641 

85,480 

1000 

756,39 

1 ,6204 

0,(4108 

81,750 

1030 

777,74 

1 ,6374 

0,41377 

78,360 

MOO 

797.78 

1 ,6527 

0,39362 

75,260 

M50 

810,65 

1.GC67 

0,37956 

72,290 

1200 

635,26 

1 ,6797 

0,37059 

69,430 

FTD-MT-24- 1924-71 
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